CONVERGENT INTERPOLATION TO CAUCHY INTEGRALS OVER 
ANALYTIC ARCS WITH JACOBI-TYPE WEIGHTS 



LAURENT BARATCHART AND MAXIM YATTSELEV 

Abstract. We design convergent multipoint Fade interpolation schemes to Cauchy 
transforms of non-vanishing complex densities with respect to Jacobi-type weights on 
analytic arcs, under mild smoothness assumptions on the density. We rely on the work 
[9] for the choice of the interpolation points, and dwell on the Riemann-Hilbert approach 
to asymptotics of orthogonal polynomials introduced in [33] in the case of a segment. 
We also elaborate on the 9-extension of the Riemann-Hilbert technique, initiated in [37] 
on the line to relax analyticity assumptions. This yields strong asymptotics for the 
denominator polynomials of the multipoint Fade interpolants, from which convergence 
follows. 



1. Introduction 

Classical Pade approximants (or interpolants) and their multipoint generalization are 
probably the oldest and simplest candidate rational-approximants to a holomorphic function 
of one complex variable. They are simply those rational functions of type'^ (m,rt) that 
interpolate the function in m+n+l points of the domain of analyticity, counting multiplicity. 
Classical Pade approximants refer to the case where interpolation takes place in a single point 
with multiplicity m -I- n -j- 1 [38] . 

Besides their everlasting number-theoretic success [43, 32, 41], they are common tools in 
modeling and numerical analysis of various fields, ranging from boundary value problems 
and convergence acceleration [16, 29, 13, 26, 20] to continuous mechanics [6, 51], quantum 
mechanics [8, 52], condensed matter physics [42], fluid mechanics [40], system and circuits 
theory [12, 31, 17], and even page ranking the Web [14]. 

In spite of this, the convergence properties of Pade or multipoint Pade approximants are 
still far from being understood. For particular classes of functions like IVIarkov functions, 
some elliptic functions, and certain entire functions such as Polya frequencies or functions 
with smooth and fast decaying Taylor coefficients, classical Pade approximants at infinity 
are known to converge, locally uniformly on the domain of analyticity [36, 50, 4, 34]. But 
when applied to more general cases they seldom accomplish the same, due to the occurrence 
of "spurious poles" that may wander about the domain of analyticity. Further distinction 
should be made here between diagonal approximants (i.e. interpolants of type (m,™)) and 
row approximants (i.e. interpolants of type (rn,n) where n is kept fixed), and we refer 
the reader to the comprehensive monograph [7] for a detailed account of many works on 
the subject. Let us simply mention that, for the case of diagonal approximants which is 
the most interesting as it treats poles and zeros on equal footing, the disproof of the Pade 
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conjecture [35] and of the Stahl conjecture [15] have only added to the picture that classical 
Pade approximants are not seen best through the spectacles of uniform convergence. 

The case of multipoint Pade approximants is somewhat different, since choosing the inter- 
polation points offers new possibilities to help convergence. However, it is not immediately 
clear how to use these additional parameters. The theory was initially developed for Markov 
functions (i.e. Cauchy transforms of positive measures compactly supported on the real line) 
showing that multipoint Pade approximants converge locally uniformly on the complement 
of the smallest segment containing the support of the defining measure, provided the in- 
terpolation points are conjugate symmetric [25]. The crux of the proof is the remarkable 
connection between rational interpolants and orthogonal polynomials: the denominator of 
the n-th diagonal multipoint Pade approximant is the rt-th orthogonal polynomial of the 
measure defining the Markov function, weighted by the inverse of the polynomial whose 
zeros are the interpolation points (this polynomial is identically 1 for classical Pade approxi- 
mants). The conjugate symmetric distribution of the interpolation points is to the effect that 
the weight is positive, so one can apply the asymptotic theory of orthogonal polynomials 
with varying weights [49] . 

When trying to generalize this approach to more general Cauchy integrals than Markov 
functions, one is led to consider non-Hermitian orthogonal polynomials with respect to 
complex-valued measures on more general arcs than segments, and for a while it was un- 
clear what could be hoped for. In the pathbreaking papers [44, 46, 47, 48], devoted to the 
convergence in capacity of classical Pade approximants to functions with branchpoints, it 
was shown that such orthogonal polynomials lend themselves to analysis when the measure 
is supported on a system of arcs of minimal logarithmic capacity linking the branchpoints, 
in the complement of which the function is single- valued. Shortly after, the same type of 
convergence was established for multipoint Pade approximants to Cauchy integrals of con- 
tinuous (quasi-every where) non-vanishing densities over arcs of minimal weighted capacity, 
provided that the interpolation points asymptotically distribute like a measure whose po- 
tential is the logarithm of the weight [27]. Such an extremal system of arcs is called a 
symmetric contour, or S-contour, and is characterized by a symmetry property of the (two- 
sided) normal derivatives of its equilibrium potential. The corresponding condition on the 
distribution of the interpolation points may be viewed as a far-reaching generalization of the 
conjugate-symmetry with respect to the real line that was required to interpolate Markov 
functions in a convergent way. 

After these works it became apparent that the appropriate class of Cauchy integrals for 
Pade approximation should consist of those taken over S'-contours, and that the interpola- 
tion points should distribute according to the weight that defines the symmetry property. 
However, it is not so easy to decide which systems of arcs are S'-contours, since finding a 
weight making the arcs of smallest weighted capacity is a nontrivial inverse problem, and in 
any case convergence in capacity is much weaker than locally uniform convergence. 

For the class of Jordan arcs, new ground was recently broken in [9] where it is shown 
that such an arc, if rectifiable and Ahlfors regular at the endpoints, is an S-contour if and 
only if it is analytic. The proof recasts the 5-property for Jordan arcs as the existence of a 
sequence of "pseudo-rational" functions, holomorphic and tending to zero off the arc, whose 
boundary values from each side of the latter remain bounded, and whose zeros remain 
at positive distance from the arc. There are in fact many such sequences that can be 
computed explicitly from an analytic parameterization of the arc. Then, translating the non- 
Hermitian orthogonality equation for the denominator into an integral equation involving 
Hankel operators and using compactness properties of the latter, the reference just quoted 
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establishes that multipoint Padc approximants to Cauchy transforms of Dini-continuous 
(essentially) non-vanishing densities with respect to the equilibrium distribution of the arc 
converge locally uniformly in its complement when the interpolation points are the zeros of 
these pseudo-rational functions. 

Still the above result remains unsatisfactory, for the hypotheses entail that the density 
with respect to arclength in the integral goes to infinity towards the endpoints of the arc, 
since so does the equilibrium distribution. In particular, ultra-smooth situations like the 
one of Cauchy integrals of smooth functions over analytic arcs are not covered. The present 
paper develops a new technique to handle any non- vanishing integrable Jacobi-type density 
under mild smoothness assumptions, thereby settling more or less the issue of convergence in 
multipoint Padc interpolation to functions defined as Cauchy integrals over analytic Jordan 
arcs. 

We dwell on the Riemann-Hilbert approach to asymptotics of orthogonal polynomials with 
analytic weights, pioneered on the line in [11, 19] and carried over to the segment in [33]. We 
also elaborate on the 9-extension thereof, initiated on the line in [37] to relax the analyticity 
requirement. This will provide us with strong (i.e. Plancherel-Rotach type) asymptotics 
for the denominator polynomials of the multipoint Pade interpolants we construct and for 
their associated functions of the second kind, from which the local uniform convergence we 
seek follows easily. It is interesting to note that the Riemann-Hilbert approach, which is 
typically a tool to obtain sharp quantitative asymptotics, is here used as a means to solve a 
qualitative question namely the convergence of the interpolants. 

The interpolation points shall be the same as in [9], namely the zeros of a sequence 
of pseudo-rational functions adapted to the arc. Such an interpolation scheme will prove 
convergent for all Cauchy integrals with sufficiently smooth density with respect to a Jacobi 
weight on the arc at the same time. This provides us with a varying weight which is not 
of power type, nor in general converging sufficiently fast to a weight of power type to take 
advantage of the results of [3], where the Riemann-Hilbert approach is adapted to non- 
Hermitian orthogonality with analytic weights on smooth S'-arcs. Instead, when "opening 
the lens" , we set up a sequence of Riemann-Hilbert problems with varying contours whose 
solutions converge to the desired one by properties of the pseudo rational functions. 

We pay special attention to keep smoothness requirements low, in order to obtain as 
general a result as the method permits. Roughly speaking, the higher the Jacobi exponents 
the smoother the density should be, see the precise assumptions (2.9). When the Jacobi 
exponents are negative, only a fraction of a derivative is needed, which compares not too 
badly with the Dini-continuity assumption in [9]. In the present setting, however, the density 
cannot vanish whereas some weak vanishing is still allowed in [9] . We are of course rewarded 
here with stronger asymptotics. 

As the varying part of our weight is analytic, the extension inside the lens with controlled 
9-estimates, introduced in [37] for power weights, needs only deal with the density defining 
the Cauchy integral we interpolate. This step is treated using either tools from real analysis, 
e.g. Muckenhoupt weights and Sobolev traces, or else classical Holder estimates for singular 
integrals, whichever yields the best results granted the Jacobi exponents. 

Since we consider analytic arcs only, it is natural to ask how general our results with 
respect to the general class of Cauchy integrals over rectifiable Jordan arcs. It turns out 
that they are as general as can be, because if the Cauchy integral of a nontrivial Jacobi 
weight can be interpolated in a convergent way with a triangular scheme of interpolation 
points that stay away from the arc, then the arc is in fact analytic [10]. 
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The paper is organized as follows. Section 2 fixes notation and defines pseudo rational 
functions as well as multipoint Pade approximants before stating the main results. In 
Section 3, the contours that will later be instrumental for the solution of the Riemann- 
Hilbert problem are introduced. Section 4 contains preliminaries on smooth extensions from 
boundary data in domains with polygonal boundaries. Section 5 is devoted to key estimates 
of certain singular integral operators that play a main role in the extension of the weight. 
Section 6 and 7 deal with the analytic Riemann-Hilbert problem, while Section 8 solves the 
d version thereof. Finally, in Section 9, we gather the material developed so far to establish 
the asymptotics and the convergence of multipoint Pade approximants stated in Section 2. 

2. Statements of Results 

Let A be a closed analytic Jordan arc with endpoints —1 and 1. That is to say, there 
exists a holomorphic univalent function S, defined in some domain D-^ D [—1, 1], such that 

A = S(hl,l]), S(±l) = ±l. 

We call S an analytic parameterization of A. We orient A from —1 to 1 and, according to 
this orientation, we distinguish the left and the right sides of A denoted by A+ and A~, 
respectively. It will be convenient to introduce two unbounded arcs, say. A; and A^, that 
respectively connect — oo to —1 and 1 to +oo, in such a manner that A;U AU A,, is a smooth 
unbounded Jordan arc that coincides with the real line in some neighborhood of infinity. 
Define on A the Jacobi weight 

(2.1) w{z)^w{a,l3;z):^{l-z)°'{l + zf, a,/3>-l, 

where we choose branches of (1 — z)" and (1 + z)^ that are holomorphic outside of A^ 
and A/, respectively, and assume value 1 at the origin. In particular, w is analytic across 
A° := A \ {±1}. Further, set 

(2.2) 115(2;) := \/ — 1, lr(z)/z^l, as z ^ 00, 
to be a holomorphic branch of the square root outside of A. Then 

(2.3) (/3(z) := z + rD(z), zeL>:=C\A, 

is holomorphic in D\ {00}, has continuous boundary values on A^, respectively, and 
satisfies 

(2.4) ip^Lp~=\ on A and </?(z)/2z— >1 as z ^ 00. 

It is immediate that is inverse of the Joukovski transformation J(z) :— (z^ + l)/2z, i.e., 
J((^(z)) = z, z G D. Moreover, maps D conformally onto an unbounded domain whose 
boundary is an analytic Jordan curve [9, Sec. 3.1] which is symmetric with respect to the 
transformation z ^ \/z. In particular, does not vanish. 

2.1. Symmetric Contours. Multipoint Pade approximants to a given function / are de- 
fined to be rational interpolants to /. In this paper we are interested in those functions / 
that can be expressed as Cauchy integrals of Jacobi-type complex densities defined on A 
(see the smoothness assumptions in (2.8) and (2.9)). In order for multipoint Pade approx- 
imants to converge to such a function, it is necessary to choose the interpolation schemes 
appropriately with respect to A. We presently characterize these schemes in terms of the 
associated monic polynomials vanishing at the interpolation points. 
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Let {vn} be a sequence polynomials such that deg(i;„) < 2?! and each Vn has no zeros on 
A. To this sequence we associate a sequence of "pseudo-rational" functions, say {r^}, given 



where the product is taken over all zeros of Vn according to their multiplicities. It is easy to 
see that each function Vn is holomorphic in D, has the same zeros as w„ counting multiplic- 
ities, and vanishes at infinity with order 2n — deg(z)„). Hence, each r„ has exactly 2n zeros 
counting multiplicities. Moreover, the unrestricted boundary values exist continuously 
from each side of A and satisfy r+r~ =1 by the first part of (2.4). 

Hereafter, the normalized counting measure of a finite set is the probability measure that 
has equal mass at each point counting multiplicities. Below, the weak* topology refers to 
the duality between complex measures and continuous functions with compact support in C. 

Definition 1. We say that a sequence of polynomials with no zeros on A belongs to 

the class S(A) if the following conditions hold: 

(1) the associated functions r„ via (2.5) satisfy \r^ \ ~ 0{1) uniformly on A and r„ — 
o(l) locally uniformly in D; 

(2) there exists a neighborhood of A that contains no zeros of r„ for all n large enough; 

(3) the normalized counting measures of zeros of r„ form a weak* convergent sequence. 

The third requirement in the definition of S(A) is purely technical and is placed only to 
simplify the forthcoming considerations since one can always proceed with subsequences as 
far as convergence is concerned. 

Regarding the nature of the class S(A), the following result was obtained in [9, Thm. 1]. 
For a closed analytic Jordan arc A, there always exist sequences {vn} belonging to S(A) and 
they can be constructed explicitly granted the parameterization S. A partial converse is also 
true. Namely, let A be a rectifiable Jordan arc with endpoint ±1 such that for x ~ ±1 and 
alH G A sufficiently close to x it holds that \At^x\ < const. |a; — , (3 > 1/2, where \At^x\ is 
the length of the subarc of A joining t and x and "const." is an absolute constant. If there 
exists a sequence of polynomials {vn} meeting the first two requirements of Definition 1, 
then A is necessarily analytic. The class S(A) is also intimately related to the so-called 
symmetry property of the contour A [44, 45, 9]. 

For our investigation we need to detail further the properties of the just defined interpo- 
lation schemes. We gather them in the following theorem. We agree that the arcs involved 
have endpoints ±1. Moreover, we say that two holomorphic functions are analytic continua- 
tions of each other if they are defined on domains that have nonempty intersection on which 
the functions coincide. 

Theorem 1. Let A be a closed analytic Jordan arc and {f„} € S(A). Then there exists a 
sequence of closed analytic Jordan arcs {A„} such that: 

(i) there exist analytic parametrizations S„ of A„ and 'E. of A such that the functions 
S„ converge to S uniformly in some neighborhood of [—1, 1] as n oo; 

(ii) for each function r„, associated to Vn via (2.5), there exists an analytic continuation 
r,*, holomorphic in Z3„ := C \ A„, such that \{rn)^\ = 1 on A„. 

Let tt)„ and ipn be defined relative to A„ as ro and (p were defined in (2.2) and (2.3) 
relative to A. Clearly, tt)„ and are analytic continuations of W and (p to D^- In fact, r* 
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is simply the function associated to w„ via (2.5) with Lp replaced by ipn- It is apparent that 
r* is nothing but the Blaschke product with respect to -D„ that has the same zeros as r„. 

2.2. Multipoint Pade Approximation. Let /i be a complex Borel measure with compact 
support. We define the Cauchy transform of /i as 

(2.6) /^(^):=y'M^, z^C\suM^A■ 

Clearly, is a holomorphic function in C \ supp(/^) that vanishes at infinity. 

Classically, diagonal (multipoint) Pade approximants to are rational functions of type 
(n, n) that interpolate at a prescribed system of 2n + 1 points. However, when the 
approximated function is of the from (2.6), it is customary to place at least one interpolation 
point at infinity so as to let the approximants vanish at infinity as well by construction. 

Definition 2. Let be given by (2.6) and be a sequence of monic polynomials, 

deg(w„) < 2n, with zeros in C\ supp(/i). The n-th diagonal Pade approximant to ff^ asso- 
ciated with {vn} is the unique rational function n„ =Pn/ln satisfying: 

• Aegpn < n, degg„ < n, and g„ ^ 0; 

• iQn{z)ff,{z) -p„(z)) /v„{z) is analytic in C\supp(/i); 

• {qn{z)ff_,{z) - pn{z)) /Vniz) ^ O {l/z'^+^) aS Z ^ GO . 

A multipoint Pade approximant always exists since the conditions for p„ and g„ amount 
to solving a system of 2n+ 1 homogeneous linear equations with 2n + 2 unknown coefficients, 
no solution of which can be such that = (we may thus assume that g„ is monic); note 
that the required interpolation at infinity is entailed by the last condition and therefore n„ 
is, in fact, of type (n ~ 1, n). 

We consider only absolutely continuous measures that are supported on A and whose 
densities are Jacobi weights (2.1) multiplied by suitably smooth non-vanishing functions. 
This leads us to define smoothness classes C™'^. 

Definition 3. Let K be an infinitely smooth closed Jordan arc or curve. We say that 
9 e C"^'''{K) if 9 is m-times continuously differentiable on K with respect to the arclength 
and its m-th derivative is uniformly Holder continuous with exponent i.e., 

<const.|ii-t2|S h,t2eK. 

When K — A, we simply write C™''' instead o/C'"'^(A). We also write C'^{K) for the 
space of infinitely differentiable functions on K . 

Together with C™'*", we also consider fractional Sobolev spaces. 

Definition 4. Let K be an infinitely smooth Jordan arc or curve. We say that 9 € 

W^^/P(if),pe(l,oo), if 

p 

\dx\\dy\ < CO. 





9{x) - 9{y) 


1 JKxK 


x-y 



When K — A, we simply write Wp instead ofWp ^^^(K). 

We shall be interested only in the case p E (2, oo) since in this range it holds that 

(2.7) Wi-i/PcCO'S ^=1-^' pe(2,oo). 

by Sobolev imbedding theorem (see Section 4.1). 



CONVERGENT INTERPOLATION AND JACOBI-TYPE WEIGHTS 



7 



In what follows, we assume that the measure /i in (2.6) is of the form 

(2.8) dfi{t) = {wh){t)dt, h{t)^e'"-*\ teA, 

where the Jacobi weight w — w{a,P; ■) and the complex function 6 are such that 

(2.9) a, /3 e (-s,s) n (-l,oo) 
with 

1-1, if 0eW^^/^ pe(2,cx3), 
2^-1, if deC°-', <;€ 

m + <;, if 6'eC™'S meN, ?e(0,l]. 

To describe the asymptotic behavior of the approximation error to functions by the 
multipoint Pade approximants, we need to introduce complex geometric means and Szego 
functions. The geometric mean of ft, = is given by 

(2.10) Gh:=exp\ [ eduj\ , duj{t) :^ , t e A. 



TlXV+{t) ■ 

The measure uj is, in some sense, natural for the considered problem as suggested by the 
forthcoming Theorem 2. Observe also that uj simply becomes the normalized arcsine distri- 
bution on A when A = [—1,1]. As J duj = 1, Gh depends only on h and is non-zero when 
is Holder continuous (see Section 5.3). Moreover, in this case the Szego function of h, given 
by 

(2.11) Sh{z) -.^ exp J ^^diuit) - ^ J edcuY zeD, 

is the unique non-vanishing holomorphic function in D that has continuous unrestricted 
boundary values on A from each side and satisfies 

(2.12) h = GhS+S^ on A and Sh{oo) = 1. 
The main result of the paper is the following theorem. 

Theorem 2. Let A be a closed analytic Jordan arc connecting ±1 and {vn} € S(A). Let 
also ff^ a Cauchy integral (2.6) with fi given by (2.8) and (2.9). Then {!!„}, the sequence 
of diagonal multipoint Pade approximants to ff^ associated with is such that 

(/^-n„)»- [2G^ + o(l)]^2^„, 

with o(l) satisfying 

(0, ™1'"'-'^'^ ), .-max{|a|,|/?|}<l, 

(0,i), s-max{|a|,|/3|} > 1, 

locally uniformly in D, where the constant const, depends on a, d^i = fidoj, i.e., /i = 
—i'Kwhvo^ , and the functions rn are associated to the polynomials w„ via (2.5) and hence 
converge to zero geometrically fast in D. 

The convergence theory of Pade approximants to Cauchy integrals is strongly interwoven 
with asymptotic behavior of underlying orthogonal polynomials that are the denominators 
of n„. In fact, it is easy to show that — n„ = Rn/Qm where Rn is a function associated 
to Qn via (2.15), and g„ satisfies the orthogonality relations of the form (2.14), (2.16) (see, 
for example, [9, Thm. 4]). Hence, Theorem 2 follows from Theorem 3 below. 



, „N 1 /, M const. 

(2.13) |o(l)|<— -, ae 
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2.3. Strong Asymptotics for non-Hermitian Orthogonal Polynomials. In this sec- 
tion we investigate the asymptotic behavior of polynomials satisfying non-Hermitian orthog- 
onality relations of the form 

(2.14) / t^qn{t)wn{t)dt = 0, j e{0,...,n-l}, 

J A 

together with the asymptotic behavior of their functions of the second kind, i.e., 

(2.15) Rn[z) -.^ R.,i{qn; z) ^ — ;, z e D, 



t — z iri 

where {w„} is the sequence of varying weights specified in (2.16). 

Theorem 3. Let {(?„}, deg((7„) < n, be a sequence of polynomials satisfying orthogonality 
relations (2.14) with weights given by 

(2.16) := , /i = e , /i„ = e 

where 6 and w = w{a,(3;-) are as in (2.9), {On} is a normal family in some neighborhood 
of A, and {vn} G S(A). Then, for all n large enough, the polynomials qn have exact degree 
n and therefore can be normalized to be monic. Under such a normalization, we have that 



(2.17) 



<Z„ - [l + o(l)]/5„ 
RntV = [l+o(l)]7„S'„ 

with o(l) satsfying (2.13) locally uniformly in D, where 

(2.18) Sn := (2/(^)"S'„^n,+ , 7„ := 2^^^"G^^„+, 

and Rn was defined in (2.15). Moreover, it holds that 

.2ln^ f 1" = [l + o{l)]/S+ + [l + o{l)]/S- 

^' ' I (i?«tr)± = [l + o(l)]7n5± 

where o(l) satisfies (2.13) locally uniformly in A° . 



The method of proof can also be used to derive the asymptotics of g„ and Rn around 
±1 as was done in [33]. However, the corresponding calculations are lengthy and do not 
impinge on the convergence of Pade interpolants proper, which is why the authors decided 
to omit them here. 

The appearance of the normal family {On} in (2.16) is not necessitated by Theorem 2 but 
is included for possible application to meromorphic approximation [53]. 

3. Proof of Theorem 1 and ^-Functions 

In this section we prove Theorem 1. The notion of function, which we introduce along 
the way, will be needed later on for the proof of Theorem 3. 

3.1. Parameterization S, functions g and g. Let {«„} G S(A) and r„ be associated to 
Vn by (2.5). As required by Definition l-(2) and (3), the normalized counting measures of 
the zeros of r„ converge weak* to a Borel measure v, supp(:/) C D. Denote by Vj^ the Green 
potential of this measure with respect to D. It was shown in the course of the proof of [9, 
Thm. 1, see (4.34)] that Definition 1-(1) yields 

Lp{z) - (p{t) 



(3.1) l/^(z) = - log 



1 - (p{z)(p{t) 



dv{t). 
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In other words, the Green kernel — log 



, where is the conformal map of D onto 



{|z| > 1} such that (j>{oo) = oo and 0'(oo) > 0, can be replaced by the one in (3.1) for this 
special measure v. 

The Green potential Vj^ is a positive harmonic function in _D \ supp(i^) whose boundary 
values vanish everywhere on A. Let Lp := {z : Vj^{z) = logp}, p > 1, be a level line of 
Vj^ in S(Z?h), the range of S. Without loss of generality we may assume that p is a regular 
value and therefore Lp is a smooth Jordan curve encompassing A = ii. Denote by O the 
domain bounded by Lp and A. Set 

(3.2) df:^]^{d^f-idyf) and df:^]^{dj + idyf). 

Since v is a. probability measure, it can be verified as in the proof of [9, Thm. 1, see (4.39) 
and after] that the function 



(3.3) m - exp |2y^ -^m] = exp |- j log f^^^^^K^) 

is well-defined in O and maps it conformally onto the annulus {z : 1 < |2;| < p} while 
$(±1) = ±1, where we take any path from 1 to z contained in O \ A. Moreover, by direct 
examination of the kernel in (3.3), we get that 

(3.4) $+ = $^=1/$" on A. 

This, in particular, yields that Jo<i> is holomorphic across A, where J{z) — (z+l/z)/2 is the 
Joukovski transformation. Consequently, the inverse (Jo <i))~^ is a holomorphic univalent 
map in some neighborhood of [—1, 1] that analytically parametrizes A. In what follows, we 
assume that S = ( J o ^)~^ . 

Based on the conformal map <&, we define two more functions, g and g as follows. Set 
L := $-i(hP,-l]). 2 := ^"H[1.P]) (see Fig. 1), and define 

5:=log$, lmi,g(z) = 0, .geH(0\L), 

^^■^^ g:=\og<^-^i, limg(z) = 0, geH(0\Z). 

2— 1 

It follows immediately from (3.4) that 

(3.6) = —g~ and Ij^ = —'g~ on A. 

Hence, g^ and are analytic in Og := (O U A) \ L and Og := (O U A) \ L, respectively. 
Moreover, it holds that ^^(A) = ^2^^^ ^ [-77^,0] and ^^(i) = g^{^l) = 0. It is also true 
that g^ and g^ are univalent in Og and Og, respectively. Indeed, suppose that g'^{zi) — 
g'^{z2), zi,Z2 € Og. Then either ^(zi) — ^{z2) and therefore zi — Z2 by conformality of <& 
or <i>(zi) = l/<i>(z2), which is possible only if $+(2:1) = <I>~(z2), i.e., if zi ^ Z2 & A. The 
case of g^ is no different. 

3.2. Jordan arcs A„, functions gn and 5„. By Definition l-(2) and upon taking p smaller 
if necessary, we may assume that functions r„ have no zeros in O. Moreover, as r„ is 
holomorphic in D \ A and has 2n zeros in D \ O, its winding number is equal to ~2n on 
any positively oriented curve homologous to Lp and contained in O. In other words, r„ 
has a continuous argument that decreases by Ann as A is encompassed once in the positive 
direction. Thus, the functions 
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are well-defined and analytic in O, where Vn is the normalized counting measure of the zeros 
of r„. Moreover, as the counting measures of zeros of r„ converge weak* to v by assumption, 
the functions <i>„ converge to $ uniformly in O, distinguishing the one-sided values on A*. 
Hence, we can define 

5„:=log$„, lini 5r„(z) = 0, .g„eH(0\L), 

(3-7) _ ^ _ _ ~ 

5„ := log$„ - TTi, lim g„(z) = 0, g„eH(0\L). 

^— 1 

By (2.4), it is straightforward to see that ^n^n = 1 on A, and therefore 
(3.8) g+ = -g,7 and g+ = -g,7 on A. 

Thus, and ^ are analytic in Og and Og, respectively. We choose domains Ol C Og and 
Oj^ C Og in such a manner that Ol 13 L, Oj; D L, and U O^;^ is simply connected and 
contains A (see Fig. 1). Then it is an easy consequence of the convergence of <&„ to $ that 
g^ and g^ converge uniformly to g^ and g^ on Ol and O-j^, respectively. 




Figure 1. The domain O bounded by A and Lp (light grey), the domain Oz, C 
O U A (dark grey), the cuts L and L (dashed arcs). 



Next, we claim that g^ and g^ are univalent for all n large enough in Ol and Oj^, 
respectively. Assume to the contrary that there exist two sequences of points {zi^„}, {z2^n} C 
Ol such that (7^(zi_„) = (7^(2:2, «)■ As Ol is compact, we can assume that Zj_„ zj € Ol, 
j — 1,2. Since g^ converges to g^ uniformly on Ol, we have that g'^{zi) — g^{z2) and 
therefore zi = Z2- Set (i„(z) := (g^(z) — (7^(zi_„))/(z — zi „). Then (i„ are analytic functions 
on Ol that converge uniformly to d{z) := (g^(z) — g'^{z\))l{z — zi). Moreover, the values 
dn{z2,n) are equal to and converge to d{zi). Thus, (5^)'(zi) = d{zi) — 0, which is 
impossible since g^ is univalent. This proves the claim as the case of is no different. 

From the above we see that each g^ maps Ol conformally onto a neighborhood of zero as 
(7^(1) = 0. Set A„ 1 to be the preimage of the intersection of this neighborhood with E2 ■= 
{C : Arg(C) = tt}. Then A„_i is an analytic arc with one endpoint being 1. Analogously, g^ 
maps Oj^ conformally into another neighborhood of zero. Thus, we can define A„ _i to be 
again the preimage of the intersection of this neighborhood with S2. Clearly, A„._i is an 
analytic arc with one endpoint being —1. Noticing that g'^ assumes negative values if and 
only if g^ assumes negative values on the common set of definition Ol H ; "^"^ derive that 
A„ :— A„_i U A„ _i is an analytic arc with endpoint ±1. 



CONVERGENT INTERPOLATION AND JACOBI-TYPE WEIGHTS 



11 



3.3. Parameterizations S„, functions 17* and g*. Now, define rD„ and cpn with respect 
to A„ like ft) and ip were defined in (2.2) and (2.3) with respect to A. Clearly, ipn is an 
analytic continuation of (p from D onto Dn = C \ A„. Further, let r* be defined by (2.5) 
with if replaced by ipn while keeping the same zeros as r„. Hence, r* and r„ are analytic 
continuations of each other defined in D and respectively. Finally, set <I>* , g*, and 
to be the analytic continuations of gn, and gn from O, O \ L, and O \ L, onto 
0„ := (O U A) \ A„, 0„ \ L, and 0„ \ L, defined in, by now, obvious manner. Hence 

3n (ffn)^ and gf^ = (g;)^ 

while (g*)^ is negative on A„. Since, in addition, (5,*)"'' = —(5^)" on A„, (g*)^ are pure 
imaginary on A„ and |<i>* | = 1 there. Therefore, (<&*)+ = (*&,*)" and (r*)+ = (r*)-. Fur- 
thermore, <&* maps 0„ onto some annular domain having the unit circle as a component of 
its boundary. Arguing as was done after (3.4), we derive that Jo 4>* is holomorphic across 
An and that S„ := (J o <&* )~^ is a holomorphic parameterization of A„, S„([— 1, 1]) = A„. 
Moreover, as the functions <i>„ converge to (f> uniformly in some annular domain encom- 
passing A, we see that the functions J o $* converge locally uniformly to J o $ in some 
neighborhood of A. Hence, the sequence of analytic parameterizations of A„ converges 
uniformly to the analytic parameterization S of A in some neighborhood of [—1,1]. This 
finishes the proof of Theorem 1 . 

4. Trace Theorems and Extensions 

As is usual in the Riemann-Hilbert approach to asymptotics of orthogonal polynomials, 
we shall need to extend the weights of orthogonality from A into subsets of the complex 
plane. As the weights are not analytic, this extension will require a special construction that 
we carry out in this section. 

4.1. Domains with Smooth Boundaries. In this section we suppose that is a bounded 
simply connected domain with boundary F which is infinitely smooth and contains A, i.e., 
A C F. 

Definition 5. Set LP(il), p e [1, 00), to be the space of all measurable functions f such that 
l/l^ is integrable over Q,. The Sobolev space Wp(f2), p G [l,cx)), is the subspace ofljP{fl) 
that comprises of functions with weak partial derivatives also in LP(f2). 

Then the following theorem takes place [30, Thm. 1.5.1.2]. 

Theorem Tl. For each f G Wp"^/^(F), p e (l,oo), there exists F £ W^(f)) such that 
F^p = /. Moreover, the extension operator can be made independent of p. Conversely, for 

every F € W^(n) it holds that F\r € wl'^^^^{T). 

Together with the Sobolev spaces Wp{il), we consider smoothness classes C™^'"(ri). 

Definition 6. By C™'^(ri), m e <; e (0,1], we denote the space of all functions on 
whose partial derivatives up to the order m are continuous on f2 and whose partial derivatives 
of order m are uniformly Holder continuous on f2 with exponent Moreover, CQ''''(f2) will 
stand for the subset o/C'"'''(f2) consisting of functions whose partial derivatives up to order 
m, including the function itself, vanish on F. Finally, C°^'{fl) will denote the space of 
functions on whose partial derivatives of any order exist and are continuous on f2. 
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It is known from Sobolev's imbedding theorem [2, Thm. 5.4] that 

' L2p/(2-p)(f]), pe [1,2), 



(4.1) wlin) C 



Hence, for / e Wp ^^^(F), p G (2, oo), the function F granted by Theorem Tl belongs to 
C0^i-2/p(n) and therefore / e C^'^-^^p{T), which is exactly what was stated in (2.7). 
Later on, we shall need the following proposition. 

Proposition 4. Let f be a continuous function on A such that /(±1) =0. If f G Wp , 
p e (2,cxd), then there exists F g Wp(ri) such that — f. Moreover, if f e 0°'*^, 
<j e (1/2, 1], i/iera i/iere exists F £ Wj(ri) /or an?/ q £ (2, jz^) smc/i that := /. 

Proof. In both cases set / = / on A and / = on r\ A. When / S Wp"^^'', it is immediate 
to check that / G Wp~^^''(r) and therefore the first claim follows from Theorem Tl. When 
/ e GO'S it holds that / G C''^'^(r) and / e Wl^^^\T) for any g e (1, j^) by an easy 
estimate (see Definition 4). Hence, the second claim of the proposition again follows from 
Theorem Tl. □ 

To state a trace theorem for classes C™''"(ri), we need to introduce the notion of a direc- 
tional derivative. Namely, let ^ be a continuous function on and / G Wp{il). With the 
slight abuse of notation, we define the derivative of / in the direction of the field ^, denoted 
by d^, as 

(4.2) d^f -.^ ^Bf + idf = Re(e)a,/ + Im(09,/, 

where V/ is the gradient of /, ^ is the vector field with values in M? corresponding to ^. 

As r is infinitely smooth, any conformal map of onto the unit disk belongs to C°° (fi) 
[39, Thm. 3.6]. Moreover, it holds that 0' 7^ in 17. Thus, we may set 

(4.3) n{z):^^ and n(z) 0(z) 

Then n,n £ C°°{^1) and n is holomorphic in fi. Moreover, for any z £ T, n{z) = n{z) /\n{z)\ 
represents the complex number corresponding to the outer unit normal to T at z. Then the 
following theorem takes place [30, Thm. 6.2.6]. 

Theorem T2(l). Let {fk}T=o that fk £ C"'-'''^T), m £ N, <i £ (0,1], k £ 

{0, . . . , to}. Then there exists F £ C"'^%n) such that (9^F)|r ^ fk, k £ {0, . . . , m}. 

Now, observe that 

fc-i 

J=l 

where the functions c^.j involve sums and products of the powers of the iterated directional 
derivatives of \(j>'\ with respect to the field n and therefore belong to C°°(ri). Set Ck,k ■= I'P'l'' 
and Ck,j = 0, j £ {k + 1, . . . , to}, k £ {1, . . . , to}. Then the matrix = [ck,j]k^j^i is such 
that det('^0) = which is non-vanishing at any z £C, and 

(a„,...,a;")T = <^49„,...,9™)T. 

Thus, to every family of functions {fk}kLoi fk G C™^'"''''(r), there corresponds another fam- 
ily, say {fk}T=o^ fk e C'"-^-^(r), such that there exists F £ C"'''^(U) satisfying (a^F)|r = A 
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and (9^F)|r — ft- Moreover, this correspondence is one-to-one and onto. Hence, Theo- 
rem T2(l) can be equivalently reformulated as follows. 

Theorem T2(2). Let {fk}T=o such that fk G C™-'^''^(r), m e N, <; e (0,1], k e 
{0, . . . ,to}. Then there exists F e C"'''Xn) such that (a^i^)|r ^ fk, k e {0, . . . ,m}. 

Finally, we define r := in. Clearly, t{z)/\t{z)\, z e F, is the complex number corre- 
sponding to the positively oriented unit tangent to F at z. Since n and r are holomorphic 
functions such that r — in, it is a simple computation to verify that dndrF — drdnF in 51. 
Then the following proposition holds. 

Proposition 5. Let f e C™^'^(A), m G N, cr € (0, 1], /^^Hil) = 0, S {0, . . . , m}. Then 
there exists F e C'^'^U) such that F\^ = / and BF e C""^'''(n). 

Proof. Set /o = / on A and /o = on F \ A. It is clear that fo £ C™'^(F). Further, set 
fk := (-i)^9^/o, fc e {1, . . .^m}. As fk G C™-^^'=(F), G {0, . . . , m}. Theorem T2(2) yields 
that there exists F G C™''^(r2) such that (a,^,F)|r = /fe. In particular, F^^ = f. 
It remains only to show that dF G C™~^'''(ri). It can be easily checked that 

(4.4) 2ndF = 9„F -I- idrF =: H in fl' , 

where f2' is an annular domain such that F C dfl' and n is non-vanishing on this domain. As 
n G C°°(n) and is zero free in n', it holds that G C"~^'''(n) if and only if 77 G €""^''(0). 
Moreover, it is immediate from the construction of F that H G C'"^^'^(ri). Thus, it is only 
necessary to verify that all the partial derivatives of H of order k, for any k G {0, . . . , to— 1}, 
vanish on F. Since partial derivatives with respect to n and r commute, and these fields are 
non- vanishing and non-coUinear in fi', it is enough to show that 

^ (a^a^^i/)|r = + ^ {d';^+'d:^-F)^^ 

for all fci + fc2 G {0, . . . , TO — 1}. The latter holds since 

(d^^d^^+'F)^^ = = (-z)'=^+i5^+'=^+Vo = = -^ (d^'+'d^F)^^ , 

by the choice of {fk}- Q 

4.2. Domains with Polygonal Boundary. The previous results also hold, with some 
modifications, for domains with polygonal boundaries. Namely, let be a domain whose 
boundary is a curvilinear polygon consisting of two pieces, say Ai and A2, such that they 
might form corners at the joints. As we do not strive for generality at this point, we assume 
that each Aj is an analytic arc connecting —1 and 1. 

The first trace theorem of this section states the following [30, Thm. 1.5.2.3]. 

Theorem T3. Given fj G Wp^^^^(Aj), j = 1,2, satisfying /i(±l) = /2(±1), there exists 
F G Wp(ri) such that F\i^. — fj, j — 1,2. The choice of F can he made in such a way that 
it depends only on fj and not on p. 

To state an analogous theorem for the classes C"''''(Aj), we again need to define normal 
fields on Aj, j = 1,2. Let F^ be an infinitely smooth Jordan curve such that A^ C Fj. 
Moreover, assume that the interior domain of Fj, say fij, contains fi, j = 1,2. Define nj for 
fij as it was done in (4.3). Composing with a self-map of the disk if necessary, we can choose 
conformal maps in (4.3) so rij does not vanish in f2. Further, set Tj := inj if lies on the 
left side of Aj and Tj otherwise. In particular, the fields Tj and rij commute, are 

infinitely smooth, non-vanishing and non-coUinear. Finally, observe that (4.4) holds with 
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n, r, and f2' replaced by Uj, Tj, f2, and the plus sign replaced by the minus sign in the 
right-hand side of (4.4) when tj = —inj. 

With all the necessary material at hand, we can state a special case of the trace theorem 
for smoothness classes on domains with polygonal boundary'^ [30, Cor. 6.2.8]. 

Theorem T4. Given {fjk}T=0' fjk e C™-'='<'"(A,), m e N, ^ e (0,1], j = 1,2, satisfying 
fjtl\±l) = 0, fci + /ca e {0, . . . , m}, there exists F e C™''^(n) such that {d^^F)\A, = fjk, 
fee {0,...,™}, i = 1,2. 

Now, as in Section 4.1, we shall make Theorems T3 and T4 suit our needs. Let A be a 
closed analytic Jordan arc and A± be two closed analytic Jordan arcs with endpoints ±1 
such that the interior domain of A U A+, say ri_|_, is simply connected and lies to the left of 
A while the interior domain of A U A_, say ri_, is again simply connected and lies to the 
right of A. Then the following proposition holds. 

Proposition 6. Let f be a continuous function on A such that /(±1) =0. /// G Wp 
p € (2, oo), then there exist F± G Wp(fl±) such that 

(4.5) F±|A = ±f and F±|a± = 0. 

Moreover, if f e C"'"- , <; G (1/2,1], then there exist F± G Wj(fi±) for any q G (2, j^) 
satisfying (4.5). 

Proof. This proposition follows from Theorem T3 in the same fashion as Proposition 4 
followed from Theorem Tl. □ 

Finally, we state the counterpart of Proposition 5 for domains with corners. 

Proposition 7. Let f G C™''=(A), m G N, c; G (0, 1], f^^\±l) = 0, A: G {0, . . . ,m}. Then 
there exist functions F± G C™''"(ri±) such that 

(4.6) f±|A=±/, J^±|A±=0, and dF± e C^-'^'in^). 

Proof. First, we consider the case of ri_|_. By setting fik '.= {—i)^f''^\ G {0, . . . we 
see that fi^ G C''^'^(A), fc G {0, . . . , m}. Moreover, after putting f2k = 0, we observe that 

/j^^^H±l)=0, fci + fc2G{0,...,m}, J = 1,2. 

Then the existence of F+ in C"'^(n+) follows from Theorem T4. The fact that dF+ G 
C™"^'''(ri+) can be shown exactly as in Proposition 5. In the case of the only difference 
is that we need to set fik :— ~i^ f'^^^ since this time the normal and tangent on A satisfy 
T = —in. □ 



5. Scalar Boundary Value Problems 
In this section we dwell on smoothness properties of certain integral operators. 



In Theorem T4 we use non-unit normal fields rij rather than fields that are unit on Aj as it was done 
in the original reference. However, we have already explained after Theorem T2(l) that these formulations 
are equivalent. 
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5.1. Integral Operators. Below we introduce contour and area integral operators and 
explain the solution of a certain 9-problem. 

Let 4> be an L^', p > 1, function on A, where L'' = L^'(A) stands for the space of functions 
with p-summable modulus on A with respect to arclength differential \dt\. The Cauchy 
integral operator on A is defined by 

(5.1) C(i>{z) := CaHz) ^ 7^ f ^^dt, zeD. 

zni J A t — z 

It is known that C4> is a holomorphic function in D with traces on A, i.e., non-tangential 
limits a.e. on A, from above and below denoted by C^(f). These traces are connected by the 
Sokhotski-Pemelj formulae [24, Sec. 1.4.2], i.e., 

(5.2) C+(/)-C"0 = and C+ cj) + cf) = Scj), a.e. on A, 
where S is the singular integral operator on A given by 

(5.3) S4>{t) 5a0(t) = - / P^dt, r e A°, 

m J^t - T 

with the integral being understood in the sense of the principal value. 

Let now f2 be a simply connected bounded domain with smooth boundary T. We define 
Cr and Sp by (5.1) and (5.3) integrating this time over F rather than A. The Sokhotski- 
Plemelj formulae (5.2) still hold for cf) G L^(r), p > 1, with the only difference that now Cr0 
is a sectionally holomorphic function and therefore Cjt </) is the trace of Cpcj) from within Q 
and Cp is the trace of Cr4> from within C \ fi. 

Concerning the smoothness oiCr4> the following is known. If g C°'^(r), € (0, 1), then 
Cr(f> e C'''''(ri) [24, Sec. 5.5.1]. In particular, this means that Cr(f> extends continuously from 

to r. Further, if is continuously differentiable on F, then Cp0 = Cr4>' [24, Sec. 4.4.4]. 
Thus, we may conclude that when G C™^'^(F), meZ+, ^ e (0, 1), then Crcf) G €'"^'■"(0). 

Let now (/> G LP(ri). The Cauchy area integral on is defined as 

(5.4) !€(/)( z) := -1^ // ^^dC AdC, zE n. 

27ri JJq Q- z 

Then it is well-known [5, Sec. 4.9] that 

(5.5) dICcj) — 4' and dICcj) — ^4'^ 

in the distributional sense, where B is the Beurling transform, i.e., 

(5.6) ^^(^^^^i/X^^^''^"'^" 

and the integral is understood in the sense of the principal value. 

The transformation K, defines a bounded operator from LP(ri) into L^P/*^^~P)(ri) for p € 
(1,2), [5, Thm. 4.3.8], and into C^^^/^lJi) for p £ (2,oo), [5, Thm. 4.3.13]. Since nothing 
prevents us from taking z outside of il, IC(j> is, in fact, defined throughout C and is clearly 
holomorphic outside of f2 and vanishes at infinity. Moreover, /C^ is continuous across F 
when p € (2, oo). The latter can be easily seen if we continue by zero to a larger domain, 
say n, and observe that this extension is in LP(ri). 

The Beurling transform B defines a bounded operator from a weighted space LfJ(C) := 
{/ '■ /^^^ G L^'(C)} into itself when the non-negative function v is an Ap-weight (Muckenhoupt 
weight), p e (l,oo) [5, Thm. 4.9.6]. Let </> e LP{n). We can suppose that (j) G Lp(C) with 
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= outside of and therefore (p/m E L^^^p{C). It holds that is an Ap-weight for 
p>2 [28, Sec. 9.1.b]. Thus, 6((/./tr) G L|'^|p(C) and therefore 

(5.7) <j>eLP{n) imphes m B{(j)/ to) E LP{n), p > 2. 

Finally, we point out that (j) E Wp(fi) can be recovered by means Cr and /C in the following 
fashion: 

(5.8) (j) ^ Cr(t) + ICd(t) a.e. in n, 
which is the Cauchy-Green formula for Sobolev functions. 

5.2. Functions of the Second Kind. Let i?„ be given by (2.15) with qn satisfying (2.14) 
and Wn defined as in (2.16). Clearly, i?„ is holomorphic in D, and it vanishes at infinity 
with order at least n + 1, i.e., Rn ~ 0(z^"^^) as z ^ cx), on account of (2.14). It is also 
clear that i?„ = 2C((7„w„). Thus, it holds by (5.2) that 

~ R~ = 2qnWn. 

Further, since qnWn/w — qnhnh/vn is Holder continuous by the conditions of Theorem 3, 
we have that 

( 0(|l-z|"), if a<0, 

(5.9) Rn=< 0(log|l-z|), if a = 0, 

[ 0(1), if a>0, 

and analogous asymptotics holds near —1. Indeed, the case a < follows from [24, Sec. 
1.8.3 and 1.8.4]. (Observe that we defined (1 - t)", t e A°, as the values on A of (1 - z)", 
where the latter is holomorphic outside of the branch cut taken along A^. However, (1 — t)" 
equivalently can be regared as the boundary values of (1 — z)" on A+, where the latter is 
holomorphic outside of the branch cut taken along A; U A. Hence, the analysis in [24, Sec. 
1.8.3] indeed applies to the present situation.) The case a — follows from [24, Sec. 1.8.1 
and 1.8.4]. Finally, the case a > holds since i?„(l) exists for such a as Wn{t)/{t — 1) is 
integrable near 1 in this situation. 

5.3. Szego Functions. Let 9 E C^'S m E Z+, ^(0, 1], and h := e" . The definition of the 
Szego function given in (2.11) can be rewritten as 

Sh = exp jmC (^^^ - i y edu;^ . 

Note that decomposition (2.12) easily follows from the Sokhotski-Plemelj formulae (5.2). 
Moreover, as the lemma in the next section shows, the traces belong to C™''" , < < 
and S^{±1) — S',~(±l). In particular, the functions'^ :— /S'f^ and := / are 
continuous on A and assume value 1 at ±1. It also follows from the Sokhotski-Plemelj 
formulae that 

(5.10) c^=exp|ir=^5 

The following facts are explained in detail in [9, Sec. 3.2 and 3.3]. First, if 6*1, 02 G C™'^, 
then Sfixh^ = ShiSh2- Second, if {On} is a normal family in some neighborhood of A then 
{Sh^} is a normal family in D. If, in addition, {9n} converges then {S'/i„} converges as well 




Here we slightly abuse the notation and use superscripts + and — as a part of the symbol for the 
function. However, in Lemma 16 we shall construct a function cjj whose traces on A will coincide with c^. 
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and the convergence is uniform on the closure of D, that is, including the boundary values 
from each side. 

Third, the uniqueness of decomposition (2.12), which was shown, for instance, in [53, eq. 
(2.7) and after], implies the following formula for the Szego function of the polynomial Vm 
deg(z;„) < 2n, with zeros in D: 

(5-11) C-^^' 

where r„ was defined in (2.5). 

Fourth, observe that it is possible to define continuous arguments of (z + l)/ip(z) and 
(z — 1) /^p{z) that vanish on the real axis in some neighborhood of infinity. Therefore it holds 
that 

(5.12) S^{z)=[2—\ f2-^j and = 2-("+«, 

where w was defined in (2.1) and the branches of the power functions are taken such that 
the positive reals are mapped into the positive reals. 

Finally, using (5.12) with w = z«(l/2, 1/2; •) we have that 

S+{t)S-{t) = 2^/1 - i2 = -2i\x>+{t) teA. 

Hence, we get that 

(5.13) S„+ = y/2m/ip and G„+ = i/2, 

where, as usual, the branch of the square root is chosen so that is positive for large 
positive reals. It will be useful for us later to note that 

(5.14) i^S^.f = iS^,)^-^SZ. = (5±^)2 J^^T^ ^ ±,sZ,, 
where we used (2.12) and (5.13). 

5.4. Smoothness of a Singular Integral Operator. In this section we show that the 
boundary values on A of the Szego function of have essentially the same Sobolev or Holder 
smoothness as 9. 

We start with the case of functions in Wp p £ (2,oo). Observe that Wg ^^"^ D 
Wp when q < p, which is immediate from Definition 4. 

Proposition 8. Let 9 £ wl^^^^, p G (2, oo). Then 

ro'^S{9/Ki+) ^±d + ni'^£, d(±l) = 0, 

where d G Wq ^^'^ for any q £ (2,p) and £ is a polynomial, deg(£) < 1. 

Proof. It follows immediately from Cauchy integral formula and the Sokhotski-Plemelj for- 
mulae (5.2) that 

(5.15) C ( ^1 = and S ( ^ | = 0. 

Hence, for any polynomial £o and t G A° it holds that 

1 = - / M^l^^^l f io{t)-£o{T) dt ^ 

^+)^^' mJ^t-TXV+{t) TTiJ^ t-T W+it) 
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where £i is a polynomial, deg{£i) < deg(^o), since is a polynomial in t. Choose 

£o to be the polynomial interpolating 6 at ±1, deg(fo) < 1- Then 



Thus, it holds by (5.2) that 



[i^J " [ V~t^/ " ^ " ^"^ ^ ^ ^^^^ ^ 
where ^2, deg(^2) < 1 will be chosen later. Set £ £i + 2^2 and 

d := 2tv+ (^C+ - ^2) -{0- io) =■■ 2di - (0 - 4). 

Clearly, to prove the proposition, we need to show that di € Wq and (ii(±l) = 0. 

Let r be any infinitely smooth curve containing A. Assume also that the inner domain 
of r, say ri, lies to the left of A, i.e., A+ is accessible from within il. Define 6e\A ■=0 — ^0 

and 6'e|r\A = 0. It is clear that 9^ £ Wp ^^^(F). Moreover, since 9^ is identically zero on 
r\ A, it holds by (5.2) that 



< const.|z2 - l|-2/p-i/2^ z^fl^ 



where from now on we agree that rD|r is the trace of tr from within f2, i.e. it is equal to ro^ 
on A. Furthermore, we can regard di as a function holomorphic in fl. Thus, by Theorem Tl, 
to show that S Wj"^''' it is enough to prove that di £ Wg(fi), q G (2,p). As di is 
holomorphic in fi, it is, in fact, sufficient to get that ddi G L'^{il). 

Now, Proposition 4 insures that there exists Q S W^(ri) such that 0|r — 9e- Observe 
that e/tv e Wl{n) for any s g [1, ^) since de/ni,de/ni G L^i^) by Holder inequality 
and Q/w^ G L^{il.) by the estimate 

ze{z) 

where we used the definition of and (4.1). Thus, Cauchy-Green formula (5.8) applied to 
Q/w implies that 

ddi ^d(e-tvic (^^^ - (0(2^ = ae - 9 (^w/c " "'^^ - ^ 

a.e. in il, where we used that dni = 0. Since dQ E L^(ri) and 1x1^2 is bounded, it is only 
necessary to show that 

belongs to L«(f]), q S (2,p), where we used (5.5). The fact that niB{dQ/ro) € LP{n) 
follows from (5.7). Now, to show that (l/tr) (/C(9e/tr) - £2) € L9(f)), q G (2,p), recall that 
ae/lr e L''(r2) for any s e [1, ^). So, as mentioned after (5.6), /C(9e/lr) e L2^/2-'*(f)) 
when 39 e (2,4], i.e., s £ (|,2); and /C(ae/rt5) G C°'i-2A(f^) when 39 e (4, 00), i.e., s can be 
chosen to lie in (2, In the first case, we get that {l/ni)IC{dQ/n>) G q £ (2,p), 
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simply by applying Holder inequality once more. This shows that ddi E L^(ri), q £ (2,p), 
when p G (2,4] with £2 = 0. In the second case, let £2 be the polynomial interpolating 
IC{de/w) at ±1, deg{£2) < 1. Then 



rv{z) 



ro 



(^) 



< const, z 



1 



(s-4)/2s 



Since 



1 2' p 



it holds that (l/tt3) (/C(ae/tt3) -4) e L'?(r2), q e (2,p), which 

shows that Sdi G L'?(ri), g e (2,p), when p e (4, 00). 

It only remains to show that di(±l) = 0. As 9^ G C'^-'^/p{T) by (2.7) and 6le(±l) = 0, 
the function C^{0(,/fo^) is either bounded near ±1 or blows up there with the order strictly 
less than 1/2 [24, Sec. 1.8.4], see also (5.9). Thus, rD+C+(6'e/rD+) vanishes at ±1 and so does 
di. □ 

We continue with the case of functions in C™''". 

Proposition 9. Let 9 e C™^^. m e Z+, e (0, 1]. Then 

tv^S{e/n>+) = ±d + n>^£, d'-''\±l)^0, k € {0, . . . ,m}, 

where £ is a polynomial, deg(^) < 2m + 1, and d G C™''' when <j G (O, ^) U (|, l), while 
d S C™'''"'^ /or arbitrarily small e > when <; = 1. 

When G (0, ^), the conclusion of the proposition follows from [21, Thm. 3]. Therefore, 
we are required"' to prove Proposition 9 only for <j G [5,1]- To this end, we shall need 
several geometrical lemmas. In all of them we assume that F is as in Proposition 8 and we 
omit superscript + for ro when dealing with the values of ro on A. By Cr we shall denote a 
constant such that |t| < Cr, r G F. Moreover, n and T2 will stand for two different points 
on F satisfying 

(5.16) |l-,i| 
Lemma 10. Let N e N or N = -1. Then 



r?|<|l-T.?| 



(5.17) 



1 



1 



< Ci|Ar|max 



tr^(Ti) n,N{T2) 
where C'l is a constant depending only on F. 

Proof. If N is an even integer, then 



1 



1 



jl-rfl 1+^/2' |l-r2|i+JV/2 



In -T2 , 



tt)^(ri) tr^(T2) 



|1 -t2|^^/2 



1 



N/2 



|l-r2|^/2|l. 

by (5.16). If N is an odd integer, then 
1 1 



Ar/2-1 

E 

3=0 



l-r 



< 



NCr\Ti-T2\ 
|1 -r2| 1+^/2 



tr^(Ti) tD~(T2) 



|1 -r2|^/2|l-r2|i/2 



Ki{t2) - tt)(ri) 



2x (A'-l)/2 



II 



,2|^/2|l_^2|l/2 



(iV-l)Cr|Ti~r2| 
|l-rf|^/2|l-r2| 



The authors were surprised not to find this case in the literature. 
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by the first estimate and (5.16). Clearly, it only remains to prove the lemma for N ~ — 1. 
It can be readily verified that it is enough to consider |ri — T2I small enough. As ro is zero 
free in fl, interior of F, an argument function, say a, is well-defined and continuous in il. 
Since tn extends holomorphically across A° and r\ A, the trace of a is uniformly continuous 
on any compact subset of F \ {±1}. Moreover, it also has one-sided limits at ±1 with the 
jumps of magnitude tt/2. Thus, there exists 6 > such that for all |ri — T2I < S it holds 
that |a(Ti) -a{T2)\ < ^. Then 

|tt(ri) +tr(T2)| > |tr(r2)| - ^-\ro{T,)\ > ^\w{t2)\ 
for \ti — T2I < (S by (5.16) for the last inequality and therefore 



|tr(ri)~m(T2)| <2Cr- 



{0{ti) + n){T2) 



< 4Cr 



|Ti 



r||l/2' 



which finishes the proof of the lemma. □ 
Lemma 11. Let g e C°^'^(r), e (i, l), and g{±l) = 0. Then (g/tv) £ C"'''-^/'^{T) and 



(5.18) 



e(n) g{T2) 



ro{Ti) ir(T2) 



< C2 min 



|l-r2|i/2'|l-r||V2 



Tl - T2 



where C'2 is a constant depending only on T. 

Proof. We start by proving (5.18). By the condition of the lemma it holds that 

\g{n) - g{T2)\ < M\n - T2\' 
\g{r)\<M\l 



(5.19) 



r2K 



Ti,T2,T e r, 



for some finite constant M. Set, for brevity, k := g/to. First, let ti — 1. Observe that 
\k{t)\ < M\1 - t^\'-^/^ < M{1 + Cr)T 



r2|l/2 



by (5.19). Thus, k{1) = by continuity and it holds that 



(5.20) 



|K(Ti)~«(r2)| <M(1 + Cr) 



In -T2r 
l-r||V2 



Clearly, an analogous bound holds when ri = —1. 

Second, let |ti — r2| > |1 — T;^|. In this case, we also have that 

(5.21) < |rf-r|| + |l-r2| <C?|ri-T2|, Q 1 + 2Cr 

Then it follows from (5.19), (5.16), and (5.21) that 

\^{n)~n{r2)\ < |^(n)|-hk(r2)|<M(|l-r^ri/2+|l 



2K-1/2 



(5.22) 



< 2M|1 - r||'^"i/2 < 2M{C^Y 



Tl -T2 



r||i/2 



Third, let \ti - T2\ < \l ~ Tf\. Then, it also holds that 



2, ^ M _2| 



|l-r2l<|l 



(5.23) 

Thus, (5.19) and Lemma 10 imply that 
Hti)-k{t2)\ < le(ri) 



w(ri) tr(T2) 



< M\l-Tt\'Ci 



T1-T2 



|g(n) - g(T2)| 

Mr2)\ 

Tl-T2\' 



1 



r2|3/2 



M 



1 



r2|l/2- 



CONVERGENT INTERPOLATION AND JACOBI-TYPE WEIGHTS 



21 



Using the conditions |ti — T2I < |1 — | and (5.23), we obtain that 



(5.24) < M(Civ/Q + l)|^_^2|i/2- 

This finishes the proof of (5.18). 

Finally, it can be readily verified that the equations leading to (5.20), (5.22), and (5.24) 
also yield that k G C^''^~^^^{T) and hence, the lemma is proved. □ 



Lemma 12. Let g be as in Lemma 11. Then 5r((?/rt') G C°'''""i/2(r) and 

1 1 

|l_^2|l/2'|l_^2|l/, 



(5.25) \Sr (I) (ri) - Sr (|) (r^) | < C, min \ „ ^ In - r^^, 



where C3 is a constant depending only on T . 

Proof. Since k := q/xv g C°'^"i/2(p)^ 5^,^ g co,?-i/2(-p) [24, Sec. 1.5.1]. To 

prove (5.25), one needs to trace the local character of the proof in [24, Sec. 1.5.1]. This is a 
tedious job but the authors felt compelled to carry it out for the reader. 
Define 

6(r) := 6rK(r) — k(t) = — / dt, t G 1 . 

m Jy- t — t 

In the light of (5.18), it is enough to show (5.25) with Srn replaced by S. Observe also that 
the integral that defines S is no longer singular as k(±1) = 0. 

Denote by V* the connected component of F n {r : |ti — r| < 2|ri — T2I} that contains ri. 
we order ti and T2 so that (5.16) holds. Then we can write 

S{r,)-S{r^) = 1/ ^^ll^^dr - 1 / ^lill^^d. 

m J J., T -T2 m Jy' t -Ti 

1 [ «(ti) - '«(t2)^^_^ f {t2 - Ti){k{t) - k{t2)) 



m Jr\r' t - n m Jr\T' {t - ti){t - T2) 

= L^ + h + Iz + Ii. 

Before we continue, observe that there exists a finite constant M such that 

(5.26) \T{t,T)\<M\t~Tl i,TGF, 

since F is a smooth Jordan curve, where |F(i,T)| is the arclength of the smallest subarc of 
F connecting t and r. 

First, let us estimate Ii. We get from (5.18) and (5.26) that 

1,1 ^ C2 . r 1 1 ^ [ \dT 

\li\ < — sup mm 



- Toll-? 



n rer- 1|1-t2|i/2' |1 - r||i/2 / j^,. \r - T2 
(5-27) < 3^ -2iT7^ / 3— ^ 



Clearly, an analogous estimate can be made for I2. 
Now, we shall estimate I3. It holds that 



(K(ri) - k(t2)) log ( — — 

n-Ti 
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where and T{, are the endpoints of F*. As |Ta — ti| — \Tb — ti\, we have that 

|ri-T2|^ 



(5.28) 1/3 1 < const. 



^|l/2' 



where const, is the product of C2 and the maximum of the argument of for all possible 

choices of and T},. 

Finally, let us estimate I4. Observe that |t — ri| < 2|r — T2I, r e F \ F*, since 

\t~ti\ < |r-r2| + |ri-T2| < |t - T2I + - n]. 
Then we get from (5.18) and the bound above that 

C2 . r 1 1 1 /■ \T2-Tl\\dT\ 



< — max mini- KrTj:^,^. ^^ttt^ | / 



2r'^ J Jr\r-' |T-ri||r-T2|i-'^ 



r\r* 



T - Tl 



T2 



I-? 



\dr\ 



(5.29) < 



C2 |t2-ti| 
^ |1-t||i/2 

C22^-^Af^-^ \t2-ti\ ds ^ C2M^-' |t2-ti|^ 



(5.30) 



Combining (5.27), (5.28), and (5.29) with (5.18), we see that (5.25) holds. □ 
Lemma 13. Let g e C°'", v G (O, 5), g{±l) = 0, be such that 

\q{ti) - g{T2)\ < C4min '^^ " ^2!"+'/', 

■'■i,''2 € F, Tl 7^ r2, where C4 is a constant depending only on F. T/ien fog G C^'^^^^/^^p), 
Fwrt/ier, fef g e C^'^(F), N eN, g'-^\±l) = 0, j S {0, . . . , N}. Then 

j {0/^''^-') e C0^-+i/2(F) 1/ u e (0, i) , 

\ (f?/m2^+l) e C0'--l/2(r) 

Proof. To verify the first claim, assume first that |ti — T2I > |1 — r^j. Since g £ C*''"(F) and 
vanishes at ±1, it holds for some finite constant M that 

(5.31) \g{r)\<M\l-Tr- 
Then we get from (5.21) and the inequality above that 

(5.32) |(tt5e)(Ti) - {Wg){r2)\ < Af (|1 - rf]^ + |1 - rfr) < Af(l + C?)|ri - T2\' . 

Assume now that |ri — T2I < |1 — Ti\. Then we get by (5.31), (5.17), and the conditions of 
the lemma that 

\{U3g){Ti)^{tvg){T2)\ < |tt3(T2)||e(ri)-e(T2)| + |(?(Ti)||w(ri)-w(T2)| 

<r n -2ii/2^ In -T2|"+^/^ , ,2it,^ |n-T2| 

^ |1 - ^2 r Ci |^_^2|l/2 + - ^1 I ,^2|l/2 

(5.33) < (C4 + CiAf)|Ti-T2r'^ + i/^ 

Equations (5.32) and (5.33) show that (tvg) G C°-''+^^^{T). 

It remains to prove (5.30). Suppose first that v £ (O, Then, by the assumptions on 
g, it holds for some finite constant M that 



(5.34) 



Ti,r2,TeF. 

£-(^)(t) < Af|l-T2|^-J+^, 
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Thus, for |ti — T2\ > |1 — |, we have from (5.34) and (5.21) that 



to 



2N-1 



(n) tr2A^-i(T2) 



< 



2N-1 



(n) 



0(^2) 



2N-1 



(t2) 



<M(l + (Cf)^)|ri-T2r+i/^ 



For \ti — T2I < |1 — r^l, we have from (5.34) that 



9(^1) 



e(n) 



2W-1 



(ri) 



^ (7) 



0(^2) 



3=0 



XV 



2N~1 



(r2) 



tt) 



< 



e(n) 



tX) 



27V-1 



(n) 



p0)(ri)(r2-ri)^' 



3=0 

< /+M|r2-Tir+i/2 
Furthermore, it holds by (5.17) and (5.34) that 

N 

g^'Hri){r2 - ny 



XV 



2N-1 



ir2) 



2^-1(^2) 

M|t2-ti|^+" 
|l-r|r-V2 



J=0 



1 



1 



ft)2JV-i(ri) ^2^-1(^2) 



E 



gO)(ri)(T2~ri)J 



ft) 



27V-1 



(n) 



N 



< (2N -1)CiMY^ 



Af+l 



< 2NCiM 



T2 - Tl 



|l-ri^r'-''+lT2-rip+i 
|l_^2|Ar+l/2 

j-i.-l/2 



|l-Tfr'J-+"|T2-rip 

|l_^2|Ar-l/2 



|t2 - rir'+i/2 < 2iV2CiAf |t2 - ri^+i 



/2 



Clearly, the case v £ (5,1) can be handled in a similar fashion. 



□ 



Proof of Proposition 9. Clearly, we need to prove the proposition only for ^ ^,1 as these 
two cases follow from the obvious inclusion C™"'""'^ C C'"''". 

Let deg(£o) < 2to + 1, be the polynomial interpolating 6* at ±1 up to and including 
the order m. Throughout the proof we assume that F is as in Proposition 8 and that 6 is 
extended to F \ A by 4. Clearly, this imphes that 6 e C™'^(F). As li := S{£o/n>+) is a 
polynomial of degree 2m, we have that 



S 



XV 



S 



4 



to 



St 



4 



to 



|A 



Assume first that to = 0. Then — £0 satisfies the conditions of g in Lemma 11 and 
therefore Lemma 12 holds with 9 — Iq in place of g. Let £2 be the polynomial interpolating 
SriiO - 4)/tt5) at ±1 and set 



(5.35) 



d := ft) 5, 



and l:=lx^i2. 



Clearly, (i(±l) = 0. Then the conclusion of the proposition follows from Lemma 13 applied 
with £> = 5r((0-4)/rt))-4. 

Assume now that to g N. Since the derivative of a singular integral is the singular integral 
of the derivative [24, Sec. 1.4.4], observe that 



5, 



(m) 



ft) 



'-In 



XV 



()ri)N 



E 

i=o 
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where Vj are polynomials. Then SrHO - io)/^) e C"'''-^^^{T) by (5.30), applied with 
g = Vj{9 — ^o)^™ N — ji and u = ^, and the fact that singular integrals preserve Holder 
smoothness [24, Sec. 1.5.1]. Thus, Sr{{0 — £o)/xo) has m continuous derivatives on F. Let 
then d and £ be defined by (5.35), where £2 is the polynomial interpolating Sr{{0 — £o)/fD) 
at ±1 up to and including the order m. Once more, since singular integral commutes with 
differentiation, we get 



- |:(:)j^(4"'-«(^)-4"-*' 



k=Q j=a 



k J \ j J ir^fe-i \ \ txi2j+i 



where Uk are polynomials and the polynomials £j_k interpolate the corresponding term in 
the parenthesis. Then 

by (5.30), applied with g = Vj+i{9 - 4) N ^ j + I, v ^ <i, and the fact that 
singular integrals preserve Holder smoothness. Thus, (5.30) applied once more, now with 
g = Sr{vj{0-£o)'^"^-''-^yn>^^+^)-ij,k, N = k, and, u = <r-l/2, yields that e C°^''(r), 
which finishes the proof of the lemma. □ 

6. RlEMANN-HlLBERT-(9 PROBLEM 

/ 0" 

In what follows, we adopt the notation 0™'^-' for the diagonal matrix I ^ 



where is a function, m is a constant, and (T3 is the Pauli matrix CT3 



1 

-1 



6.1. Initial Riemann-Hilbert Problem. Let be a 2 x 2 matrix function and Wn be 
given by (2.16). Consider the following Riemann-Hilbert problem for ?V (RHP-^^): 

(a) is analytic in C \ A and lim ^^(z)z^"'^^ — y, where is the identity matrix; 

1 2?ii„ 



(b) has continuous traces from each side of A° , , and ^ — 

(c) 3^ has the following behavior near z — 1: 

„ / 1 |1 \ 
Oi ^ ' ' , If «<0, 



1 



^/l logjl ^1 ] if „^o, as Dbz^I; 
o(^ \ \y if a > 0, 

(d) has the same behavior when I? 9 z ^ — 1 as in (c) only with a replaced by /? 
and 1 — z replaced by 1 + z. 

The connection between RHP-^^ and polynomials orthogonal with respect to w„ was first 
realized by Fokas, Its, and Kitaev [22, 23] and lies in the following. 
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Lemma 14. Let qn be a polynomial satisfying orthogonality relations (2.14) and i?„ be the 
corresponding function of the second kind given by (2.15). Further, let q*^_i be a polynomial 
satisfying 

ql_^{t)wn{t)dt = 0, j e {0, . . . , n - 2}, 



A 



and = Rn-iisin-i'^ ') function of the second kind for q*i_i. If a solution of 

RHP-^^ exists then it is unique. Moreover, in this case deg(g„) — n, i?*_j(z) = 0(z~") as 
z — > oo, and the solution of RHP-^^ is given by 



(6.1) ^ 



where m„ is a constant such that m„i?*_2^(z) = z "[1 + o(l)] near infinity. Conversely, if 
deg(g„) = ri and i?*_]^(z) = 0(z~"') as z oo, then '3/ defined in (6.1) solves RHP-^^ . 

Proof. As only the smoothness properties of the function Wn were used in [33, Lem. 2.3], 
this lemma translates without change to the case of a general closed analytic arc and yields 
the uniqueness of the solution of RHP-^^ whenever the latter exists. 

Suppose now that the solution, say ^ = [^fcl^^^i, of RHP-^ exists. Then = z"+ 
lower order terms by the normalization in RHP-^^(a). Moreover, by RHP-^^(b), has no 
jump on A and hence is holomorphic in the whole complex plane. Thus, is necessarily a 
polynomial of degree n by Liouville's theorem. Further, since = 0{z~'^~^) and satisfies 
RHP-^(b), it holds that = 2C{3^uWn). From the latter, we easily deduce that 
satisfies orthogonality relations (2.14). Applying the same arguments to the second row of 
iV, we obtain that ^^21 = 9n-i and ^2 = mnRn-i with m„ well-defined. 

Conversely, let deg(g,i) = n and i?*_]^(z) = 0(z^") as z 00. Then it can be easily 
checked by direct examination of RHP-^^(a)-(d), using the material in Section 5.2, that 
given by (6.1), solves RHP-^. □ 

6.2. Renormalized Riemann-Hilbert Problem. Throughout, unless specified other- 
wise, we follow the convention = y^[zfexp{iArg(z)/2}, Arg(2;) e (— 7r,7r]. Set 



(6-2) Sn :— \J G(„„//ift„)/2 and En :— en(p'' S^y^/hh^)- 

Then _E„ has continuous boundary values on each side of A that satisfy 

E+ e:^ = — — = — 

" " 2hh„ '2wn 
due to (2.4), (2.12), and (2.16). Further, put 

(6.3) c+:=S'+/5,7, c+ 5,+ /S"^^, c" := l/c+, and c" := l/c+. 

Then we get on account of (2.12), (5.11), and the multiplicativity property of the Szego 
functions that 

(a A) i^-^^fi^y _ _ „ „^+ and ^ - (r c c)- 

where we slightly abuse the notation by writing (r„c„c)* instead of r^c^c^. Since any 
Szego function assumes value one at infinity and ip(z)/2z — > 1 as z — > cx), it holds that 
£'„(z)/[e„(2z)"] ^ 1 as z 00. Then it is a quick computation to check that 

f 1 2Wn \ fj^+\-<^3 _ ( (r„C„c)+ W 



^^"^ ^ 1 ; ^^"^ V {rnCnCY 
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and 

lim (2"e„)'"^^£;-'"^'(z) = J. 

z — >oo 

Suppose now that RHP-^^ is solvable and '¥ is the solution. Define 
(6.5) (2"e„)'"^^£:,7'"^ 
Then solves the following Riemann-Hilbert problem (RHP-^): 

(a) 9" is analytic in D and ^(oo) = ^\ 

(b) !J has continuous traces, 3'±, on A° and = 3"- 

(c) S/' has the following behavior near z — 1: 







w 



sr ^ < 



o 
o 
o 



log 1 1 — z 
log 1 1 — z 



1 1 
1 1 



if a < 0, 
if a = 0, as D 3 z ^ \\ 
if a > 0, 

— 1 as in (c) only with a replaced by /3 



(d) has the same behavior when D 3 z 
and 1 — z replaced by 1 + z. 
Trivially, the following lemma holds. 

Lemma 15. RHP-J?' is solvable if and only if RHP-^^ is solvable. When solutions of 
HHP -.9 and lUIP-iV exist, they are unique and connected by (6.5). 

6.3. Opening the Lenses, Contours T,ext, S„, and S™''. As is standard in the Riemann- 
Hilbert approach, the second transformation of RHP-^^ is based on the following factoriza- 
tion of the jump matrix in RHP-^(b): 

(r„c„c)+ w \ f 1 '^V'^ ^ ^ 

(r„c„c)" J \ (r„c„c)"/w 1 / V / V (?'„c„c) + /w 1 

where we took into account that r^r^ = 1 on A. This factorization leads us to consider a 
new Riemann-Hilbert problem with three jumps on a lens-shaped contour S„ (see Fig. 2). 
However, to proceed with such a decomposition, we need to extend and to the complex 




Figure 2. The contour E„ := A„+ U A U A„_ C E(Ds) (solid lines). The 
extension contour 'Eext A+ U A U A_ (dashed lines and A). 



plane. We shall do it in such a manner that the extended functions, denoted by c and c^, 
are analytic outside of a fixed lens Sext (see Fig. 2). We postpone this task until the next 
section and describe here the construction of the lenses T,ext and S„. 
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We start from Y^ext- When A = [—1,1], fix a; > and set to be the subarc of the 
circle {z : \z — ix\ — |a; + 1|} that lies in the upper half plane. Clearly, A_|_ joins —1 and 1 
and can be made as uniformly close to [—1, 1] as we want by taking x sufficiently large. We 
set A_ to be the reflection of A+ across the real axis. We also denote by 51+ and ri„ the 
upper and the lower parts of the lens Sea;*, i-e., 51+ (resp. il_) is a domain bounded by A+ 
(resp. A_) and A. When A is a general closed analytic arc parametrized by S, set Y^ext to 
be the image under S of the corresponding lens for [—1,1] (the latter can always be made 
small enough to lie in i^s). 

We continue by constructing the lens S, which will we the limit lens for the sequence 
{!]„}. Let and Og C S(Dh) be defined in (3.5). Then 1 G Og, g^{l) = 0, and 
is conformal in Og. Set Us '■— {z : \z — 1\ < S}. Choose (5o > to be so small that 
Us C Og and g'^{Us) is convex for any 6 < 6o- We require the same conditions to be fulfilled 
by So and Us {z : |z + 1| < S} with respect to g^ and Og also defined in (3.5). Fix 
S < Sq. Let Jordan arcs Kj, j = 1,3, be the preimages of Si :— {( : Arg(C) — 27r/3} and 
S3 :— {( : Arg{() — — 27r/3} under 5^ in Us- Let also Kj, j — 1,3, be the preimages of Si 
and S3 under g^ in Us- Set KiU K2 U K3, where K2 is the image under S of the line 

segment that joins 'E.~^{Ki n Us). Set also 'E.~^{K^ n Us), and := U K4 U Ki, where 
K4 is the image under S of the line segment that joins S^^(if3 n Us) and 'E.^^{Ki n Us). 
Then A± are Jordan arcs that with endpoints ±1. We define S :— AUK+UK- (see Fig. 3). 

Let gn and g„ be defined by (3.7). Assume that 6 is small enough that Us C Ol and 
Us C Oj^. Then we construct the lens S„ := A U A„+ U A„_ exactly as we constructed S 
only with g, g, and S replaced by and S„, where we also employ the notation A„± 

for the upper and lower lips of the lens. It can be easily seen that the arcs A„+ and A 




Figure 3. Contours E (dashed lines) and E™'' (solid lines). Neighborhoods Us 
and Us (disks around ±1). 

intersect only at ±1 for all n large enough since A„ approach A in a uniform manner by 
Theorem 1 and A„± and A„ form angles 7r/3 at 1 and —1 by construction. 

Finally, it will be useful for us later to define one more system of contours, say S™''. The 
lens S™'' is obtained from S„ simply by replacing A by A„ (see Fig. 3). We also require 
the lens Se^t to be contained within each lens S„ (see Fig. 2). 

6.4. Extension with Controlled d Derivative. Without loss of generality we may as- 
sume that S„ C -Dh and all the functions hn are holomorphic in D^.. By the very definition 
of we have that 
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Thus, there is a natural holomorphic extension of each given by 

(6.6) -.^ Gh,^Slh-^ in D^\A. 
Concerning the extension of c, we can prove the following. 

Lemma 16. Let 6 e Wp"^^^, p £ (2,oo), or 9 e C™'S m e Z+, <r G (0, 1], m + > i . 
T/ien f/iere exists a function c, continuous m C \ A and up to A*, satisfying 

c\A±—c^, c = exp{lTi^} in C\{fl+Ufl-), and dc~cf, 

where I is a polynomial, deg(^) < 2to + 1, / G U'{^±) when 9 e wl ^'"^ , f G L«(r2±), 
g e (2, when 9 e C"^S and / e C"~^'''""(n^) when 9 e C™'S m € N, and e e (0, <;). 

Proof. This lemma is a straightforward consequence of (5.10), inclusion C"'"^ C ^'"^ for 
q £ {2, 1^) ' Propositions 8 and 9 combined with Propositions 6 and 7. □ 

6.5. Formulation of Riemann-Hilbert-i9 Problem. In this section we reformulate RHP- 
3^ as & Riemann-Hilbert-9 problem. In what follows, we understand under c and c„ the 
extensions obtained in Section 6.4 above. Suppose that RHP-=^ is solvable and 5^ is the 
solution. We define a matrix function on C \ S„ as follows: 

(6.7) ^ := I ^ ( TTnlnc/w 1 ) ' ™ 

[ 5^, outside the lens E„, 

where the upper part, ri„+, (resp. lower part, ri„_) of the lens S„ is a domain bounded by 
A„_|_ (resp. A„_) and A. This new matrix function is no longer analytic in general in the 
whole domain D since c may not be analytic inside the extension lens '2'ext- Recall, however, 
that by the very construction, c coincides with a holomorphic function c = exp{lT)i'} outside 
the lens Ti^xt- To capture the non-analytic character of , we introduce the following matrix 
function that will represent the deviation from analyticity: 

\ . „ 

TrnCndc/w y ' 

Q ^ ' outside the lens Tiext- 

Then 5^ solves the following Riemann-Hilbert-9 problem (RH9P-J^): 

(a) .5^ is a continuous matrix function in C \ E„ and ^(00) = J^; 

(b) =5^ has continuous boundary values, y±, on E° := S„ \ {±1} and 



(6.8) := 



- ^- ( -1% ) 

(c) For a < 0, has the following behavior near z — 1: 

y{z)^o(^ I jll^jl ) , as C\E„9z^l. 
For a = 0, has the following behavior near z — 1: 

^ ' \ log|l-z| log|l-0| I ^ 
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For a > 0, has the following behavior near z — 1: 

1 outside the lens E„, 
1 inside the lens S„; 
-1 as in (c) only with a replaced by 

p and L — z replaced by i + z; 
(e) =5^ deviates from an analytic matrix function according to dS^ = ^Wq. 

Then the following lemma holds. 

Lemma 17. RH9P-,5^ is solvable if and only if RHP-^ is solvable. When solutions of 
URdP-J^ and RHP-^ exist, they are unique and connected by (6.7). 

Proof. By construction, the solution of RHP-i7 yields a solution of KUdP-S^. Conversely, 
let S^* be a solution of RH9P-^. It is easy to check using the Leibnitz's rule that 9,5^* is 
equal to the zero matrix outside of S„, where 3^* is obtained from o5^* by inverting (6.7). 
Thus, is an analytic matrix function in C \ E„ with continuous boundary values on each 
side of 5]° . Moreover, it can be readily verified that has no jumps on A°_|_ and therefore 
is, in fact, analytic in D. It is aslo obvious that it equals to the identity matrix at infinity 
and has a jump on A described by RHP-5^(b). Thus, J7* complies with RHP-5^(a)-(b). 

Now, if a,/? < then it follows from RH9P-^(c)-(d) and (6.7) that has the same 
behavior near endpoints ±1 as Therefore, =5^* solves RHP-,^ in this case. When 

either a or /3 is nonnegative, it is no longer immediate that the first column of i^* has 
the behavior near ±1 required by RHP-^(c)-(d). This difficulty was resolved in [33, Lem. 
4.1] by considering ^*^~^, where ^ is the unique solution of RHP-5^. However, in the 
present case it is not clear that such a solution exists (see Lemma 14). Thus, we are bound 
to consider the first column of by itself. 

Denote by and the 11- and 21-entries of . Then £^{1 and £^2*1 analytic 
functions in D with the following behavior near 1: 

a < 
a = 0, 

a > and z is inside the lens, 
a > and z is outside the lens, 

for j — 1,2. The behavior near —1 is identical only with a replaced by /? and 1 — z replaced 
by 1 + z. Moreover, each ^'^^ solves the following scalar boundary value problem: 

(6.10) 0+ = 0-(r„c„c)+ on A, q^eB.{D). 

Now, recall that r^r^ = 1 on A and r„ has 2n zeros in D that lie away from the lens S„. 
Hence, the argument of increases by 27m when A is traversed from —1 to 1. Moreover, 
for and each a branch of the argument can be taken continuous and vanishing at 
±1 (it is the imaginary part of tv~^S{6/ni^), which is continuous and vanishing at ±1 by 
Propositions 8 and 9). Define g :— log(r„c„c) + , ^?(— 1) ~ 0. This normalization is possible 
since r+(— 1) = 1 as is a product of 2n factors each of which is equal to —1 at —1. 
Furthermore, this normalization necessarily yields that ^^(1) = 2Trni and that the so-called 
canonical solution of the problem (6.10) is given by [24, Sec. 43.1] 

0,(z):=(z-l)-"exp{C(e;z)}, z e D. 



y{z) = 




as z 



(d) =5^ has the same behavior when C \ I]„ 9 z 



(6.9) 



0(1), if 

0(log|l-z|), if 

0(|l-z|-"), if 

0(1), if 
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Recall that (pc is bounded above and below in the vicinities of 1 and —1, has a zero of order 
n at infinity, and otherwise is non-vanishing. Hence, the functions (pj := j — li2, 

are analytic in C\{±1}. Moreover, according to (6.9), the singularities of these functions at 
1 and —1 cannot be essential, they are either removable or polar. In fact, since ipjiz) = 0(1) 
or 4>j{z) = 0(log |1 ± z\) when z approaches 1 or —1 outside of the lens, (pj can have only 
removable singularities at these points. Hence, (j>j{z) = 0{1) and subsequently ,5^*]^ = 0(1) 
near 1 and —1. Thus, ^* satisfies RHP-J?'(c)-(d) for all a and /3, which means that ^* 
is the solution of RHP- 5^. Therefore, indeed, the problems RHP- 5^ and RH^P-^ are 
equivalent. □ 



Elaborating on the path developed in [37], we put KRdP-S^ aside for a while and consider 
an analytic approximation of this problem. In other words, we seek the solution of the 
following Riemann-Hilbert problem (RHP-^) : 

(a) ^ is a holomorphic matrix function in C \ S„ and ^(oo) = 

(b) £/ has continuous traces, £/±, on S° that satisfy the same relations as in RH9P- 



(c) the behavior of near 1 is described by RH(9P-^(c); 

(d) the behavior of near —1 is described by RH9P-.5^(d). 

Before we proceed, observe that the function c coincides on A„-|- with the analytic function 
c := exp{rD£}, where ^ is a polynomial, by construction. Hence, we can assume that the 
jump matrix in RHP-j2/(b) is expressed in terms of c rather than c. 

7.1. Modified RHP-^. The problem above almost falls into the scope of the classical 
approach to asymptotics of orthogonal polynomials. We say "almost" because it is not 
generally true that the functions r„ can be written as the 2n-th power of a single function, 
even up to a normal family as is the case in [3, Thm. 2] . This will explain why we constructed 
another lens, S™'', in Section 6.3. 

Consider the following Riemann-Hilbert problem (RHP-^): 

(a) ^ is a holomorphic matrix function in C \ S™'' and 3^{oo) = 

(b) has continuous traces, S^±, on (E™'^)" that satisfy 



(c) the behavior of SS near 1 is described by RHP-j2/(c) with respect to the lens S™''; 

(d) the behavior of SS near —1 is described by RHP-^(d), again, with respect to E™''. 
In fact, this new problem is equivalent to RHP-jz/. 

Lemma 18. The problems RHP-jz/ and RHP-^ are equivalent. 

Proof. Suppose that RHP-i^ is solvable and ^ is a solution. As before, let ri„+ (resp. ri„_) 
be the upper (resp. lower) part of the lens S„. Analogously define fi™^ and set 



7. Analytic Approximation of RHdP-y 



^(b); 




on a°+ua: 



(7.1) 
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Observe that fln± H fl"^:^ is a finite, possibly empty, union of Jordan domains by analyticity 
of A and A„. It is a routine exercise to verify that £/* comphes with RHP-^(a) and (b). 
Moreover, within iln+ H i^^^ and ri„_ n we have that for a < 0, £/* has the following 
behavior near 2=1: 

as z ^ 1; for a = 0, ^* has the same behavior near z = 1 as ^ since the latter is multiplied 
by a bounded matrix near 1; for a > 0, £/* has the following behavior near z — I: 

^ > 1 J^\\l-z\-" J 1 

Hence, £/* has exactly the behavior near 1 required by KKP-jz/ (c). In the same fashion one 
can check that jz/* satisfies RHP-j2/(d) and therefore it is, in fact, a solution of RHP-^. 
Clearly, the arguments above could be reversed and hence each solution of RHP-jz/ yields a 
solution of RHP-^. □ 

Let us now alleviate the notation by making a slightly abuse of it. Throughout this sec- 
tion, we shall understand under ip, r„, g„, g„, c„, c, to, Sh„, 5*1^+, and their holomorphic 
continuations that are analytic outside of A„ rather than A. Note that outside the bounded 
set with boundary A U A„ these continued functions coincide with the original ones. More- 
over, their values considered within the interior domain of A„ U A can be obtained through 
analytic continuation of the original functions across A. 

7.2. Auxiliary Riemann-Hilbert Problems. In this subsection we define the necessary 
objects to solve RHP-^. This material essentially appeared in [33] for the case A = [— 1, 1]. 

7.2.1. Parametrix away from the endpoints. As r„ converges to zero geometrically fast away 
from An, the jump matrix in RHP-^(b) is close to the identity on A°_|_ and A°_. Thus, 
the main term of the asymptotics for in D„ = C \ A„ is determined by the following 
Riemann-Hilbert problem (RHP-^): 

(a) c/K is a holomorphic matrix function in Z3„ and .yK(oo) = J^; 

w 
-1/w 

It can be easily checked using (5.14) that a solution of RHP-,yK when w = 1 is given by 

s-i ^i^s^+r' \ 

Then a solution of RHP-^/K for arbitrary w is given by 
(7.3) ^ = {G^r^'.yKiG^Slr''^/^. 

7.2.2. Auxiliary parametrix near the endpoints. The following construction was introduced 
in [33, Thm. 6.3]. Let la and be the modified Bessel functions and Ha^ and Ha"^ be 
the Hankel functions [1, Ch. 9]. Set to be the following sectionally holomorphic matrix 
function: 

nc)-nca)-( ^ 



(b) ^ has continuous traces, on A° and .yf^ = ./(^ 



(7.2) = 
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for |Arg(C)| <27r/3; 

/ iffW (2(-C)^/^) iifi^) (2(-C)V2) ^ . 

■ 1^ (i?^)) (2(-C)i/2) ttC^/^ (h!^^) (2(-C)1/^) ' 

for 27r/3 < Arg(C) < tt; 

_,cV2(/fC3)y(2(-c)V2) ,^1/2 (^wy (2(_c)V2) 1^^' 

for — TT < Arg(C) < — 27r/3, where Arg(^) € (— 7r,7r] is the principal determination of the 
argument of (. Let further Si, S2, and S3 be the rays {( : Arg(() = 27r/3}, {( : Arg((^) = 
tt}, and {( : Arg(^) = — 27r/3}, respectively, oriented from infinity to zero. Using known 
properties of /„, K^, Ha \ Ha \ and their derivatives, it can be checked that ^' is the 
solution of the following Riemann-Hilbert problem RHP-\[': 

(a) ^' is a holomorphic matrix function in C \ (Si U S2 U S3); 

(b) ^ has continuous traces, '^±, on S^, j e {1,2,3}, and 

*+ = ( ? ) on SJ, 

*+ = ( I ) oil ^2, 

*+ = ( ^-Li 5 ) on S§; 

(c) ^ has the following behavior near 00: 

r ^ 72 I i + o(c-^/2) i + o(c-i/2) 



2C^/V3 



uniformly in C \ (Si U S2 U S3); 
(d) For q; < 0, \E' has the following behavior near 0: 

|^|a/2 |^|a/2 

icr/' ici 

For a = 0, ^' has the following behavior near 0: 



* = 0( /2 ) asC-0; 



V log ICI log Id J 

For a > 0, ^' has the following behavior near 0: 

■ |^|a/2 |^|-a/2 ■ 




|^|a/2 |^|-a/2 j as C ^ in | Arg(C) | < 27r/3, 
|J|-a/2 |J|-a/2 j as C - in 27r/3 < |Arg(C)| < tt. 

$ := (T3*(-; /3)(T3, 

then this matrix function satisfies RHP-\I/ with a replaced by [3 and reversed orientation for 
S,-,ie {1,2,3}. 
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7.2.3. Parametrix near 1. As shown in [33], the main term of the asymptotics of the solution 
of RHP-^ near the endpoints of A„ is described by a solution of a special Riemann-Hilbert 
problem for which ^' will be instrumental. Let Us be as in the construction of I]„ (see 
Section 6.3). Below we describe the solution of the following Riemann-Hilbert problem 
(RHP-^): 

(a) is a holomorphic matrix function in Us \ SJ^'^; 

(b) has continuous boundary values, 3^±, on Us H (S™'')" and 

^+ = -/ on C/5n(A°+UA°_), 

= ■^-( ) 

(c) ^^-1 = + 0{l/n) uniformly on dUs] 

(d) For a < 0, ^ has the following behavior near z ~\: 

For a = 0, ^ has the following behavior near z ~ 1: 

V as us\^-'^^^l; 

\ iog|i-z| iog|i-2| y ' 

For a > 0, ^ has the following behavior near z — I: 




1 outside the lens E;"'', 



as z ^ 1 inside the lens S 



md 



To present a solution of RHP-^, we need to introduce more notation. Denote by Ug' 
and U^ the subsets of Us that are mapped by into the upper and lower half planes, 
respectively. Without loss of generality we may assume that functions 9n are holomorphic 
in Us and the branch cut of w in Us coincides with the preimage of the positive reals under 
g^. In particular, we have that w is analytic in U^ and U^ and therefore across A°_|_. Set 

^^^^^ ^ exp{i(M^ _ ^ 

\/(^h„Sh„{z) 

where we use the same branch of {z + 1)^^^ as in definition of w and a branch of (z — 1)"/^ 
analytic in Us\ An and positive for positive reals large enough. Then 

'"^'w/CnC, in Ul, 

-^°'"w/cnC, in U^, 

by the definition of c and on account of (6.6). Moreover, it readily follows from (7.4) and 
(6.3) that 

(7.5) A+A-^w on A°. 

1 /2 

Observe that gn is a holomorphic function on Us \ \i such that 

(7.6) (^1/2)+ ^ ,(^1/2)- on A„ 
by (3.8). Then the following lemma holds. 



(7.4) Al 
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Lemma 19. A solution of RHP-^ is given by 

Proof. Except for some technical differences, the proof is analogous to the considerations in 
[33, eqn. (6.27) and after]. First, we must show that fS° is holomorphic in Us- This is clearly 
true in Us \ A„. It is also clear that (o has continuous boundary values on each side of A° . 
Since 

- ,^.«r^(i, 7)(™r./^((.;/f)" 

where we used RHP-^(b), (7.5), and (7.6), S' is holomorphic across A°. Thus, it remains 
to show that has no singularity at 1. For this observe that 

since g'^^ has a simple zero at 1. Furthermore, by the very definition it holds that 

Finally, (^„/5^)(z) 2^'^"+^'^/^ as z ^ 1 by (5.12). Hence, the entries of S' can have at 
most square-root singularity at 1, which is impossible since is analytic in Us \ {1}, and 
therefore co is analytic in the whole disk Us- 

The analyticity of implies that the jumps of ^ are those of ^' (n^g^/Aj A^'^^rZ^^'^ . 
Clearly, the latter has jumps on E™'' n Us by the very definition of g^ and ^. Moreover, 
it is a routine exercise, using RHP-vE'(b) and (7.4), to verify that these jumps are described 
exactly by RHP-,^(b). It is also clear that RHP-^(a) is satisfied. Further, we get directly 
from RHP-^'(c) that the behavior of on dUs can be described by 

4 J [^ng,.) [, + oil/n) l + 0(l/n)j'^" ' 

where the property 0{l/n) holds uniformly on dUs- Hence, using that the diagonal matrices 
A~'^^ and r^^^^ commute, we get that 

= J^A^^ {J + (1/n)) A-^^'jV^^ = ^ + O (1/n) 

since the moduli of all the entries of ^ A'^ are uniformly bounded above and away from 
zero on dUs- Thus, RHP-^(c) holds. Finally, RHP-^(d) follows immediately from RHP- 
^'(d) upon recalling that |g^(z)| = 0(|1 - z\) and |A„(z)| - |1 - as z ^ 1. □ 
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7.2.4. Parametrix near —1. In this section we describe the solution of the Rieniann-Hilbert 
problem that plays the same role with respect to —1 as RHP-^ did for 1. Below we describe 
the solution of the following Riemann-Hilbert problem (RHP-^): 

(a) ^ is a holomorphic matrix function in Us \ S™''; 

(b) ^ has continuous boundary values, 3^±, on Us H and 

^+ = -/ on c75n(A°+UA°_), 

- ■^-( ) 

(c) ^jV-'^ = J + 0(l/n) uniformly on dUs\ 

(d) For /3 < 0, ^ has the following behavior near z = — 1: 

^^^(l lltzH' ^ lJs\Y.T'3z^-\- 
For /? = 0, ^ has the following behavior near z = — 1: 

^=of!°sj;t1 !°'!ItlV as c/Asr3^--i; 
\^ log |i + z| log |i + z| y ' ^ " 

For /9 > 0, ^ has the following behavior near z = — 1: 

O^l 1^' asz-^— 1 outside the lens S™'', 

/ 11 + zh^^ 1 \ 

^ ( jl _l_ ^l"'^ 1 / ' as z — 1 inside the lens S™ . 
This problem is solved exactly in the same manner as KHP-^. Thus, we set 

where the branch of (1 — z)"/^ is the same as the corresponding one in w, and (-1 - z)'^/2 
is holomorphic in Us \ A„ and positive for negative real large enough. As in (7.4), we have 
that 

e^'^'-w/cnC, in 
e~-P^^yj/cnC, in U^ , 

where Uf have the same meaning as in the previous section. However, here one needs to 
be cautious since g„ reverses the orientation on E2, i-e. S2 is now oriented from zero to 
infinity, and therefore is mapped into {z : Im(z) < 0} and U^ into {z : Im(z) > 0}. 
Again, it can be checked that 

AtK^w on A°. 

1/2 ~ 

The following lemma can be proven exactly as Lemma 19 using that (— l)"r„ = e""". 
Lemma 20. The solution of RHP-,f3^ is given by 

Finally, we are prepared to solve RHP-^/. 
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7.3. Solution of KHP-^/. Denote by T,'^ the reduced system of contours obtained from 
E™*^ by removing A„ and A„± n (Ug U Us) and adding dUs U dUg. For this new system of 
contours we consider the following Riemann-Hilbert problem (RHP-^): 

(a) ^ is a holomorphic matrix function in C \ 'Elf and ^(oo) = 

(b) the traces of are continuous on except for the branching points of 'Ejf, 
where they have definite limits from each sector and along each branch of EJj''. 
Moreover, ^± satisfy 

on dUs, 
on dUs 

.y/-^ on Y7^f\{dUs\^dUs). 
Then the following lemma takes place. 

Lemma 21. The solution of RHP-^ exists for all n large enough and satisfies 

(7.7) ^=J^ + 0{l/n), 
where 0{l/n) holds uniformly in C. Moreover, det(^) — 1. 

Proof By RHP-^(c) and RHP-^(c), we have that RHP-^(b) can be written as 

(7.8) ^+^^-{^ + 0{l/n)) 

uniformly on dUs U dUs- Further, since the jump of ^ on Y^^f \ {dUs U dUs) is analytic, 
it allows us to deform the problem RHP-^ to a fixed contour, say S'"'', obtained from S 
like T,'^'^ was obtained from E™'' (the solutions exist, are simultaneously unique, and can 
be easily expressed through each other as in (7.1)). Moreover, by the properties of r„, the 
jump of ^ on E'"'' \ {dUg U dUg) is geometrically uniformly close to J^. Hence, (7.8) holds 
uniformly on S'"'*. Thus, by [18, Cor. 7.108], RHP-^ is solvable for all n large enough 
and ^± converge to zero on E'''^ in L^-sense as fast as 1/n. The latter yields (7.7) locally 
uniformly in C \ E''^. To show that (7.7) holds at z G E'"'*, deform E'"'^ to a new contour 
that avoids z (by making 6 smaller or chosing different arcs to connect Ug and Ug). As the 
jump in RHP-^ is given by analytic matrix functions, one can state an equivalent problem 
on this new contour, the solution to which is an analytic continuation of However, now 
we have that (7.7) holds locally around z. Compactness of E'"'^ finishes the proof of (7.7). 

Finally, as 3^, and ^ have determinants equal to 1 throughout C [33, Rem. 7.1], 
det ^ is an analytic function in C \ Y.^'^ that is equal to 1 at infinity, has equal boundary 
values on each side of E'''* \ {branching points}, and is bounded near the branching points. 
Thus, det(^) = 1. □ 

Finally, we provide the solution of RHP-^. 

Lemma 22. The solution of RHP-^/ exists for all n large enough and is given by (7.1) 
with 

{.^.yy, m C\(UgUUgUY„), 
in Ug, 
m Ug, 

where Si is the solution of Moreover, det jz/ = 1. 



1 

rnCnC/w 1 
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Proof. It can be easily checked from the definition of 3^, and o/K, that given by 
(7.9), is the solution of RHP-^. As det^ = det(^) = det(^) = 1, it holds that det{£/) = 
det(^) = 1 in C, which finishes the proof of the lemma. □ 

8. 9-Problem 

In the previous section we completed the first step in solving KHdP-S^. That is we solved 
RHP-ja/, the problem with the same conditions as in RH9P-^ except for the deviation from 
analyticity, which was entirely ignored. In this section, we solve a complementary problem, 
namely, we show that a solution of a certain 9-problem for matrix functions exists. Set 

(8.1) W ^ .siW^si-^ , 

where was defined in (6.8) and ^ is the solution of RHP-jz/. In what follows, we seek 
the solution of the following 9-problem (9P-f^): 

(a) ^ is a continuous matrix function in C and i^(cxi) — J'\ 

(b) deviate from an analytic matrix function according to d'Sl = , where the 
equality is understood in the sense of distributions. 

Then the following lemma holds. 

Lemma 23. // (2.9) is fulfilled, then is solvable for all n large enough. The solution 

& is unique and satisfies 

(8.2) &=y + o{l), 
where o(l) satisfies (2.13). 

Proof. We start by examining the summability and smoothness of the entries of W. As 
£/ is the solution of RHP-^, it is an analytic matrix function in ri± and its behavior 
near ±1 is given by RH9P-c^(c)-(d). Since detsi/ = 1 in C, the behavior of jz/^^ near 
±1 is also governed by the matrices in RH9P-o?^(c)-(d) with the elements on the main 
diagonal interchanged. Observe also that |c„c| is uniformly bounded in fl±. Then a simple 
computation combined with Lemma 16 yields that 

(8.3) <const.|r„/,,^/|, ;,fc=l,2, 

where / comes from the decomposition of Be in Lemma 16 and fa.piz) := log^ |1 — if 
V := max{|Q;|, \(3\} = and fa,i3{z) ■= |1 — z'^\^' otherwise. 

Let s be as in (2.9) and denote by fl the union fi+ U fi- U A°. When 6 e Wp"^^^, 
p e (2, oo), it holds that / £ LP{il) by Lemma 16. Then we get from the Holder inequality 
that 



(8.4) f^^iif eL^in), qe 



,11 I , " u < 1, 

' 1 + 1! — S ' ' ' 

(2, oo] , s — V > 1, 



since fa,p G L«(r2), q G (^, f ). When 9 £ C°'S G (^,1], (8.4) can be obtained from 



1-? 

1-5 f ^ r^m—\fi 



the inclusion C"-'^ C Wg ^^'^ for q e (2, When 9 S C™^S m e N, <^ e (0, 1], we have 



that / G Cq (ri), e e (0,1;), by Lemma 16. This, in particular, implies that 

|(/a,/3/)(z)|<C0nst.|z^-lr-l-''-^ 

which consequently yields (8.4). 



'^Recall that by Lemma 16, / and all its partial derivatives up to and including the order m — 1 have 
well-defined vanishing boundary values on A and therefore / is indeed C"^^^''' throughout £7. 
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Suppose now that 9P-^ is solvable and ^ is a solution. Let F be a smooth arc en- 
compassing n. Convolve '3 with a family of moUifiers so that the result is smooth and 
converges in the Sobolev sense to 3) . Then by applying the Cauchy-Green formula (5.8) to 
this convolution and taking the limit, we get that 

'3 = Cr{3) + /Cr( W) = ^ + IC^3 

since 'W has compact support il, i.e., ^ is analytic outside of 17, and !^{oo) = , where 
A^^(') = IC{-W). Hence, every solution of dP-3 is a solution of the following integral 
equation 



(8.5) 



where 2 is the identity operator. As explained in Section 5.1, 1 — /C^ is a bounded operator 
from L°°(C) into itself that maps continuous functions into continuous functions preserving 
their value at infinity. Conversely, if is a solution of (8.5) in L°°(C^^^) then ^ is, in 

2 X 2 — 

fact, continuous in C , analytic outside of fi, ^(oo) = J^, and 83 = QiW by (5.5) in 
the distributional sense. Thus, dV-S) is equivalent to uniquely solving (8.5) in (L°°(C))^^^ 
because 2l — K-^Ql is holomorphic in C and is identity at infinity. 
We claim that 



(8.6) 



< 



Co 



a e 



2 > 



(0^ 



w > 1, 



where || • || is the norm of K.^/ as an operator from L°°(C) into itself and the constant Co 
depends on a. Assuming this claim to be true, we get that (X — /C^)~^ exists as a Neumann 
series and 



,1-11^)^1 

which finishes the proof of the lemma granted the validity of (8.6). Thus, it only remains to 
prove estimate (8.6). To this end, observe that (8.3), (8.4), and the Holder inequality imply 
(8.7) 



^#'11 < C\ max 



< C2 max 

1,S1 zGO 













z — ■ 


q,n 





[1,2), 



v<l, 
V > I, 



where Ci , C2 are constants depending on q and 
it holds that 



\q,n 



is the usual norm on L'^{il). Thus, 




by Holder inequality, where C3 is a constant depending on q. 

In another connection, let ip be the conformal map of D onto D, 1/1(00) — 0, tp'{oo) > 0, 
and let &„ be a Blaschke product with respect to D that has the same zeros as r„ counting 
multiplicities, i.e., 

6„(.)= n ^^^^-^^'^ 

r„(e)=0 



■ i^ie)ip{z) ' 



ze D. 



Then by the maximum modulus principle for analytic functions and Definition 1-(1), we 
have that 



(8.9) 



|r„(z)| < max|r±(t)||6„(z)| < Ci\b„{z)\, 



zeD, 
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where C4 is mdependent of n and z. Denote by Lp, p G (0,1), the level line of ip, i.e. 
Lp := {z E D : \4'{z) \ = p}- Due to Definition l-(2), there exist 1 > po > Pi > such that 
f2 is contained within the bounded domain with boundary Lp„ , say f2pp , and all the zeros of 
bn are contained within the unbounded domain with boundary Lp^ . Then 

(8.10) \K\U,n < ||&„||,,o,„ = (bn o i;-') 

where Ap^^i := {z : po < \z\ < 1}, -0^^ is the inverse ofip, and the index q is as in (8.8). As 
ip~^ is a conformal map of D onto D, it holds that |('0~"'^)'| < C5 in Ap^^i. Set 

6:(z):=(6„oV.-i)(z)= n f^' 

1 — zc 

Then by (8.8), (8.9), and (8.10) and after, we get that 

(8.11) ||/C^|| <C5||6:il,,A,„,,, 

where the constant C5 depends on q. Observe now that for \z\ — p, p E (po, 1), it holds that 
(8.12) 

^ n ^ "p{-(i -.)E ^} ^ »p{-.„<i^tei 

since |e*| < pi by the definition of pi. Clearly, (8.11) and (8.12) yield that 

< C5 (^j^J exp !^-2nq^^—^^^—^^ pdpdt^ < C^n-^'^ 
which is exactly (8.6) with a=\/q. □ 

9. Solution of RHP-^ and Proof of Theorem 3 

In this last section, we gather the material from Sections 6-8 to prove Theorem 3. It is 
an immediate consequence of Lemmas 14, 15, and 17, combined with Lemmas 18, 22, and 
23 that the following result holds. 

Lemma 24. // (2.9) is fulfilled, then the solution of 1^SP-W uniquely exists for all n large 
enough and can he expressed by reversing the transformations ^ ^ ^ using (6.5) and 
(6.7) with — '3is:^ , where ^ is the solution of RHP-i/ and ^ is the solution of dP-^. 

9.1. Asymptotics away from A, formula (2.17). We claim that (2.17) holds locally 
uniformly in D. Clearly, for any given closed set in D, it can be easily arranged that this set 
lies exterior to the lens S^'', and therefore to the lenses S„ and Sea;*. Thus, the asymptotic 
behavior of '3^ on this closed is given by 

due to Lemma 24, where e„ and £"„ were defined in (6.2), is the solution of RHP-^ given 
by Lemma 22, and ,y\^ is the solution of RHP-^ given by (7.3). Moreover, we have that 

(9.1) ^^f?-([l + 0(l)]^fc),,^l_2' 

where o(l) satisfies (2.13) locally uniformly in C \ {±1}, including on (A^)°, on account of 
(7.7) and (8.2). Thus, it holds that 

r ^Tn - [1 + o(l)]^n£;n/(2"e„) = ^"/(2«5^„+5(;,^„/„„)) 

\ ^ri2 = [1 + o(l)]^i2/(£;„2"e„) =2zG„G(,,^„/„„)5^5(,,^„/„„)/(^"+i^„+) 
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by (6.2) and (7.3), where o(l) satisfies (2.13) locally uniformly in D. Recall now that the 
entries of ,yV are, in fact, continued Szego functions defined with respect to A„. However, 
we have already mentioned that they coincide with S„+ and outside of a set exterior to 
A„ U A. Thus, the equations above indeed hold true. Hence, asymptotic formulae (2.17) 
follow from (6.1), (2.18), and (5.13). 

9.2. Asymptotics in the Bulk, Formula (2.19). To derive asymptotic behavior of qn 
and Rn on A \ {±1}, we need to consider what happens within the lens T^ext and outside 
the disks Us and Us- We shall consider the asymptotics of W from within ri_|_, the upper 
part of the lens Yiext, the behavior of '3^ in fi- can be deduced in a similar fashion. 

Recall that A„ either coincides with A or intersects it at finite number of points, as both 
arcs are images of [—1, 1] under holomorphic maps. Set 

a; := a n a;* := a n a^- , 

where and ri„_ are the upper and lower parts of the lens S'j"*. Then, it holds that 

r on a;, 

(9.2) H ^- f ' 

y \ — 1/w 

by (7.1), where with a slight abuse of notation we denote by ^± the values of in ri„± 
and on A^. Then it holds on A° by RHP-c/K(b) and on account of Lemma 24, (9.2), and 
(7.9) that 

on a: 



on a;*, 




where, again, under JV^ we understand the values of jy in ^n± and on A^ and is the 
analytic continuation of jV that satisfies RHP-^/T, only with a jump across A. Clearly, ,y\^ 
is defined by (7.2) and (7.3), where and are the Szego functions of ro^ and w with 
respect to A and not the continued functions that actually appear in (7.2) and (7.3). Thus, 
we deduce from Lemma 24 and (9.1) that 

- p"'.-)-°-(|i + ''(i)i'C),,,^,,(„..J,./„ \ )(Etr 

- (2-^j-(ii+»(i)i..«),,^,,,(^f;„ ,;^.) 

where o(l) satisfies (2.13) locally uniformly on A° and we used (6.4) to obtain the second 
equality. Therefore, it holds that 

r = [1 + o(l)] Jf+i?+/(2"e„) + [1 + 0(<5„,)]^+i?-/(2"e„«;) 

\ r+ =[l + o(l)]^+/(i?+2"6„) 

with o(l) satisfying (2.13) locally uniformly on A°. As in the end of the previous section, 
we deduce (2.19) from (6.1), the formulae 



and by noticing that 



2"e„ E+2''en 



CONVERGENT INTERPOLATION AND JACOBI-TYPE WEIGHTS 



41 



on A° by (2.12) and (5.14). 
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